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ABSTRACT 

OH 231.8+4.2, a bipolar outflow around a Mira-type variable star, displays a unique molecular richness amongst circumstellar en¬ 
velopes (CSEs) around O-rich AGB and post-AGB stars. We report line observations of the HCO^ and H'^CO^ molecular ions and 
the first detection of SO^, N 2 H^, and (tentatively) HsO^ in this source. SO^ and HsO^ have not been detected before in CSEs around 
evolved stars. These data have been obtained as part of a full mm-wave and far-IR spectral line survey carried out with the IRAM 
30 m radio telescope and with Herschel/HIFI. Except for HsO^, all the molecular ions detected in this work display emission lines 
with broad profiles (EWHM~ 50-90 km s“'), which indicates that these ions are abundant in the fast bipolar outflow of OH 231.8. The 
narrow profile (EWHM~14km s“') and high critical densities (>10® cm“^) of the HsO^ transitions observed are consistent with this 
ion arising from denser, inner (and presumably warmer) layers of the fossil remnant of the slow AGB CSE at the core of the nebula. 
Erom rotational diagram analysis, we deduce excitation temperatures of rex~ 10-20 K for all ions except for H^O^, which is most 
consistent with rex~100 K. Although uncertain, the higher excitation temperature suspected for H^O^ is similar to that recently found 
for H 2 O and a few other molecules, which selectively trace a previously unidentified, warm nebular component. The column densities 
of the molecular ions reported here are in the range Vtot~[l-8]x 10*^ cm“^, leading to beam-averaged fractional abundances relative to 
H 2 of X(HCO+)~10-*, X(H'3CO+)~2xlO-^ X(SO+)«4xlO-^ X(N2H+):«2xlO-^ and A:(H3O+)~7 x 10-‘^ cm^^. We have performed 
chemical kinetics models to investigate the formation of these ions in OH 231.8 as the result of standard gas phase reactions initiated 
by cosmic-ray and UV-photon ionization. The model predicts that HCO^, SO^, and HsO^ can form with abundances comparable to 
the observed average values in the external layers of the slow central core (at ~[3-8]xl0'® cm); HsO^ would also form quite abun¬ 
dantly in regions closer to the center (X(H3O^)~10“^ at ~10*® cm). Eor N 2 H^, the model abundance is lower than the observed value 
by more than two orders of magnitude. The model fails to reproduce the abundance enrichment of HCO^, SO^, and N 2 H^ in the 
lobes, which is directly inferred from the broad emission profiles of these ions. Also, in disagreement with the narrow H 30 ^ spectra, 
the model predicts that this ion should form in relatively large, detectable amounts (a; 10“*) in the external layers of the slow central 
core and in the high-velocity lobes. Some of the model-data discrepancies are reduced, but not suppressed, by lowering the water 
content and enhancing the elemental nitrogen abundance in the envelope. The remarkable chemistry of OH 231.8 probably reflects 
the molecular regeneration process within its envelope after the passage of fast shocks that accelerated and dissociated molecules in 
the AGB wind ~800 yr ago. 
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1. Introduction 

OF! 231.8+4.2 (hereafter, OH 231.8) is a well-studied bipolar 
nebula around an OH/IR source'. Although its evolutionary 
stage is not clear because of its many unusual properties, it is be¬ 
lieved to be a planetary nebula (PN) precursor probably caught 
in a short-lived transitional phase. The obscured central star, 
named QXPup, is classified as M9-10III and has a Mira-like 
variability consistent with an evolved asymptotic giant branch 
(AGB) star (Cohen 1981; Feast et al. 1983; Kastner et al. 1992; 
Sanchez Contreras et al. 2004). The evolution of this object may 
have been complex since it has a binary companion star (of type 


' OH/IR objects are infrared-bright evolved stellar objects with a dense 
envelope showing prominent OH maser emission. 


AO V), which has been indirectly identified from the analysis of 
the stellar spectrum reflected by the nebular dust (Cohen et al. 
1985; Sanchez Contreras et al. 2004). The system, located at 
~1500 pc (Choi et al. 2012), has a total luminosity of ~10 "'Lo 
and its systemic velocity relative to the Local Standard of Rest 
(LSR) is y^5'^~34km s '. OH231.8 is likely a member of the 
open cluster M 46 with a progenitor mass of ~3 Mq (Jura & Mor¬ 
ris 1985). 

OH 231.8 is surrounded by a massive (~ 1 Mq) and pre¬ 
dominantly cold (~10-40K) molecular envelope well character¬ 
ized by mm-wavelength emission from CO and other molecules 
(Fig. 1). The molecular gas is located in a very elongated and 
clumpy structure with two major components: (/) a central core 
(referred to as clump 13) with an angular diameter of ~6-8", a 
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total mass of ~0.64 Mq , and low expansion velocity of yexp~ 6 - 
35 km s'; and (if) a highly collimated ^6"x57" bipolar outflow, 
with a mass of ~Q.3Mq and expansion velocities that increase 
linearly with the distance from the center, reaching values of up 
to ~200 and 430 km s * at the tip of the north and south lobe, 
respectively. The temperature in the lobes is notably low, ~10- 
20K (Sanchez Contreras et al. 1997; Alcolea et al. 2001). 

The molecular envelope of OH 231.8 is markedly different 
from the slow, roughly round expanding CSEs of most AGB 
stars; its pronounced axial symmetry, large expansion velocities, 
and the presence of shocks are common in objects that have left 
the AGB phase and are evolving to the PN stage, the so-called 
pre-PNe (Neri et al. 1998; Bujarrabal et al. 2001; Castro-Carrizo 
et al. 2010; Sanchez Contreras & Sahai 2012). It is believed 
that the nebula of OH 231.8 was created as the result of a huge 
mass-loss event that occurred during the late-AGB evolution of 
the primary at a rate of MsslO^'^M© yr“*. With a total linear 
momentum of ~27 Mokm s“*, the bipolar flow is interpreted as 
the result of a sudden axial acceleration of the envelope. The lin¬ 
ear distance-velocity relation observed in the CO-outflow (with a 
projected velocity gradient of Vu~6.5km s“' arcsec“*) suggests 
that such an acceleration took place ~ 800 yr ago in less than 
~150yr. The low-velocity, low-latitude central core of the out¬ 
flow is thought to be the fossil remnant of the AGB star’s CSE. 

It is probable that the acceleration of the bipolar lobes re¬ 
sulted from the violent collision of underlying jets (emanating 
from the stellar companion) and the slowly expanding AGB en¬ 
velope (Sanchez Contreras et al. 2000; Alcolea et al. 2001; Bu¬ 
jarrabal et al. 2002; Sanchez Contreras et al. 2004); this is one 
plausible scenario that has been proposed to explain the shap¬ 
ing and acceleration of bipolar pre-PNe and PNe (e.g., Sahai 
& Trauger 1998; Balick & Prank 2002). Recently, Sabin et al. 
(2014) have found indications of a well-organized magnetic field 
parallel to the major axis of the CO-outflow of OH 231.8 that 
could point to a magnetic outflow launching mechanism. Alter¬ 
natively, as noticed by these authors, the magnetic field could 
have been dragged by the fast outflow, which may have been 
driven by a different mechanism. Por example, the underlying 
jets could have been launched by the main-sequence compan¬ 
ion powered by mass accretion from the mass-losing AGB star 
through an PU Ori type outburst (Sanchez Contreras et al. 2004) 
or a much more violent ($100 days) Intermediate-luminosity 
Optical Transient (IPOT) like event (Soker & Kashi 2012). 

OH 231.8 has a remarkably rich and unusual chemistry 
amongst CSEs around O-rich low-to-intermediate-mass evolved 
stars. In addition to the typical oxygen-rich content, with 
molecules such as H 2 O, OH or SiO (Bowers & Morris 1984; 
Morris et al. 1987; Zijlstra et al. 2001; Sanchez Contreras et 
al. 2002; Desmurs et al. 2007), OH 231.8 displays strong lines 
of many different molecular species, including many containing 
carbon. High-angular-resolution mapping of the HCO'^ (7=1-0) 
emission (Pig. 1) indicates that this ion is comparatively more 
abundant in the fast lobes that in the slow central core (Sanchez 
Contreras et al. 2000). Single-dish maps of the SiO (7=5^) 
emission show that the abundance of this molecule could also 
be enhanced in the lobes (Sanchez Contreras et al. 1997). The 
spectrum of OH 231.8 is unusually rich in lines from S- and N- 
bearing molecules. Por example, it was the first O-rich CSE 
in which H 2 S, NS, CS, and OCS were detected, and we have 
recently reported the first detection of HNCO, HNCS, HC 3 N, 
and NO (Velilla Prieto et al. 2014). Some of these S- and N- 
compounds are present in the envelope at levels not expected in 
O-rich CSEs around low-to-intermediate mass stars. In the case 
of OH 231.8, it has been proposed that extra Si and S are released 


into the gas phase from the sputtering of dust grains by shocks. 
Shocks might also initiate (endothermic) reactions that trigger 
the N and S chemistry and could also be additional suppliers of 
free atoms and ions (Morris et al. 1987; Sanchez Contreras et al. 
2000; Velilla Prieto et al. 2014). 

We have recently completed a sensitive molecular line sur¬ 
vey of this object in the mm-wavelength and far-IR range with 
the IRAM 30m telescope and Herschel (Velilla Prieto et al. 2014; 
Sanchez Contreras et al. 2014, full survey data to be published by 
Velilla Prieto et al., Sanchez Contreras et al., in prep.). We have 
detected hundreds of molecular transitions, discovered >30 new 
species (including isotopologues), and extended the rotational 
ladders for many others. This has led to very detailed informa¬ 
tion on the global physicochemical structure of this envelope. 
Prom a preliminary analysis of the survey data (Sanchez Contr¬ 
eras et al. 2014), we find two main temperature components: a) 
a predominantly cold (~10-40K) component mainly traced by 
CO (but also by, e.g., HCN, HNC, HNCO, and HCO+) and b) 
a warm (ss 100 K) molecular component that is selectively traced 
by certain molecules, such as H 2 O, amongst a few others (e.g., 
CS, H 2 S, and SiS). 

In this paper, we report the first detection of the molecular 
ions SO^, N 2 H+, and (tentatively) HsO^ in OH 231.8 as part of 
our surveys. We also present and analyse several lines of HCO'^ 
and H'^CO^, which are known to be present in OH 231.8 from 
earlier works (Morris et al. 1987; Sanchez Contreras et al. 1997, 
2000). Molecular ions are believed to be significant contrib¬ 
utors to the molecule formation process in circumstellar envi¬ 
ronments, which are active sites of molecular synthesis. In the 
interstellar medium (including molecular clouds, star-forming 
regions, PDRs, etc.) more than 30 molecular ions have been 
detected, however, in the envelopes around low-to-intermediate 
mass evolved stars detections still remain scarce and limited to 
HCO^ in most cases (see, e.g., Cernicharo & Bachiller 2011). 
In only a few objects in the post-AGB or PN phase, other pos¬ 
itive ions have been detected, namely, CO^, N 2 H'^, CH'^, OH'^, 
and (tentatively) HCS^. Except for OH^, so far these ions have 
been exclusively identified in C-rich objects, which show in gen¬ 
eral a richer chemistry than their O-rich analogues. In particular, 
toward the young PN NGC7027, where HCO+, H'3C0+, CO+, 
N 2 H^, CH^, and, tentatively, HCS'^, are observed (Cernicharo et 
al. 1997; Hasegawa & Kwok 2001; Zhang et al. 2008, and refer¬ 
ences therein) and the pre-PN CRP618, where HCO'^ and N 2 H'^ 
are detected (Bujarrabal et al. 1988; Pardo et al. 2007). Recently, 
the ion OH'^ has been detected with Herschel in five PNe (Etx- 
aluze et al. 2014; Aleman et al. 2014). To our knowledge, in cir¬ 
cumstellar envelopes negatively charged molecular anions have 
been found to date only in the C-rich AGB star IRCH-10216 (e.g., 
McCarthy et al. 2006; Cernicharo et al. 2007; Agundez et al. 
2010 a). 

2. Observations 

These data are part of a spectral line survey in the mm/far-IR 
wavelength range carried out with the IRAM 30 m radio tele¬ 
scope (~79-355GHz; Sanchez Contreras et al. 2011; Velilla 
Prieto et al. 2013) and the Herschel space observatory (~479- 
1244GHz; Sanchez Contreras et al. 2014). The spectra are pre¬ 
sented in units of antenna-temperature (T*^, which can be con¬ 
verted to a main-beam temperature (Tmb) scale via 
where rj^s is the frequency dependent ratio between the main- 
beam efficiency (rfrab) and the forward efficiency ( 771 ) of the tele¬ 
scope (values are provided in § 2.1 and § 2.2 for IRAM 30 m and 
Herschel, respectively). 
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In both surveys, we have observed one single posi¬ 
tion toward the center of OH 231.8 (with J2000 coordinates 
R.A.7H2“16?830; Dec. =-14°42'52''10). Considering the half 
power beam width (HPBW) of the IRAM 30 m and Herschel 
telescopes at the observed frequencies (HPBW 30 m ~l"-29" and 
HPBWhso ~18"-41"), this pointing fully covers the slow cen¬ 
tral core at all frequencies and the fast bipolar lobes up to some 
extent depending on the frequency (Fig. 1). 



R.A. offset (arcsec) 

Fig. 1. Distribution of the HCO^(l-O) (gray scale and solid con¬ 
tours) and '^CO(l-O) (dotted contours) velocity-integrated emission 
in OH 231.8 from previous interferometric observations; adapted from 
Fig. 2 of Sanchez Contreras et al. (2000). The smallest and largest 
HPBW of the IRAM 30 m (long dashed) and Herschel (dotted) tele¬ 
scopes in these observations are represented by the circumferences. We 
labeled the different regions/clumps in the molecular outflow as 11 to 17, 
as in Alcolea et al. (2001). Clump I2d corresponds to the region where 
the HCO^ emission peaks. The Vlsr range (in km s“') of each clump 
is 11) [-80:-30], 12) [-30 :-h 10], I2d) [-20 :-h 60], 13) [-Hl0:-t55], 14) 
[-t55:-H80], 15) [-t80:-Hl50], I6a) [-tl50:-t205], I6b) [-h205:-h 230], and 
17) [-1-230:-1-285]; see Table2 in Alcolea et al. (2001), and also Sanchez 
Contreras et al. (1997, 2000), for more detail on the physical properties 
of the clumps. 


2.1. Observations with IRAM 30 m/EMIR 

Our mm/sub mm-wave length survey was performed with the 
IRAM 30 m radio telescope (Pico Veleta, Granada, Spain) using 
the new generation heterodyne Eight Mixer Receiver (EMIR). 
Spectra were taken in several observational campaigns be¬ 
tween years 2009 and 2013. We covered the whole accessi¬ 
ble frequency range ~79-355GHz, using the four EMIR bands 
E090=3mm, E150=2mm, E230= 1 mm, and E330=0.8mm 


(Carter et al. 2012). The EMIR was operated in single sideband 
(SSB) mode for band E150, and in dual sideband (2SB) mode 
for bands E090, E230 and E330. In all cases, the image side¬ 
band was rejected with an average sideband rejection of ~14 dB. 
Each receiver band was connected to different backends; here 
we report data observed with the WILMA autocorrelator, which 
provides a spectral resolution of 2 MHz, and the East Eourier 
Transform Spectrometer (ETS) in its 195 kHz spectral resolu¬ 
tion mode. Observations were done in wobbler switching mode 
with a wobbler throw of 120". Pointing and focus were checked 
regularly (every ~1.5 and ~4h, respectively) on strong nearby 
sources. On-source integration times per tuning step were typ¬ 
ically ~ 1 h. Calibration scans on the standard two load system 
were taken every ~ 18 min. Errors in the absolute flux calibra¬ 
tion are expected to be :$25%. The parameters of the beam of 
the IRAM 30 m telescope used in this work are described to a 
good accuracy as a function of the frequency (v) by 


( HPBW(")= 

\ /7eff = 0.9exp-(™)2 


( 1 ) 


according to measurement updates performed in August 2013^. 

We reduced the data using CLASS^ following the stan¬ 
dard procedure, which includes killing bad channels, subtracting 
baseline, and averaging individual, good quality scans to pro¬ 
duce final spectra. 


2.2. Observations with Herschel/HIFI 

Observations were carried out with the Herschel Space Obser¬ 
vatory (Pilbratt et al. 2010) and its HIEEs wideband spectrom¬ 
eter (WBS de Graauw et al. 2010) in several runs in November 
2011 and April-May 2012. Our survey, executed in spectral scan 
mode, covers the frequency range ~479-1244GHz (bands la- 
5a) with an average spectral resolution of Av=l.l MHz. Obser¬ 
vations were performed in the dual beam switching (DBS) mode 
with a 3' chop throw. The two orthogonal receivers of HIEI (hor¬ 
izontal H, and vertical V) were used simultaneously. 

We reprocessed raw HerschellRWl data from level 0 to level 
2 running the hifiPipeline task of HIPE^ (versions 9.0-11.0; 
Ott 2010). Afterward, level 2 spectra were saved to PITS format 
and then imported by CLASS where standard data reduction rou¬ 
tines were applied: blanking bad-quality data, removing spurs, 
fitting and subtracting baselines, stitching of the spectrometer 
subbands, and combining individual spectra. We averaged the 
spectra from both V and H polarisations (with equal weights), 
reducing the noise in the final product. 

HIEI is a double-sideband (DSB) heterodyne instrument and, 
therefore, every spectrum contains the upper and lower sideband 
data folded together. We performed the deconvolution of the 
DSB data into single-sideband (SSB) format in CLASS using 
the task deconv. We assumed a side-band gain ratio of one. 
We performed the DSB-convolution using all bands simultane¬ 
ously, which takes advantage of redundant observations in the 
frequency overlap regions between bands and optimizes the rms 

^ http: //WWW. iram.es/IRAMES/mainWiki/EmirforAstronoiners 
^ CLASS is a worldwide software to process, reduce, and analyze 
heterodyne line observations maintained by the Institut de Radioas- 
tronomie Millimetrique (IRAM) and distributed with the GILDAS soft¬ 
ware, see http: //www. iram.fr/IRAMFR/GILDAS. 

HIPE is a joint development by the Herschel Science Ground Seg¬ 
ment Consortium, consisting of ESA, the NASA Herschel Science Cen¬ 
ter, and the HlFl, PACS, and SPIRE consortia. 
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of the final SSB spectra. The absolute flux calibration uncer¬ 
tainty of the full, deconvolved SSB spectrum is estimated to be 
~10% (HIFI Observers’ Manual, version 2.4, § 5.7). Here we 
adopt a more conservative value of ~20%. 

We adopted the values of HPBW and ? 7 eft most recently up¬ 
dated and reported by Mueller et al. (2014)^, which supersede 
previous estimates in Roelfsema et al. (2012). HIFI has 14 
mixers (7 frequency bands in 2 polarizations), each of which 
has a frequency-dependent antenna pattern (beam). Detailed 
models of the beam have been obtained independently for each 
mixer and separately for polarizations H and V (i.e., a total of 
7x2x2=28 beam models). The HPBW of HIFI for one spot 
frequency per mixer and the adopted values for the Ruze-like 
scaling of the main-beam efficiency (r/mb) within the frequency 
range of each mixer are provided in Table 3 and Table 2, respec¬ 
tively, of Mueller et al. (2014). For convenience, we reproduce 
the values of these beam-model parameters in Table A. 1 of our 
Appendix A. Taking these values into account, the HPBW and 
r/eff can be described by 


HPBW(") = 


_ 20950 


v[GHz] 

?7eff = ^ ^ exp [- 


'.2 


)1 


( 2 ) 


with the forward efficiency, rji, being 0.96 at all frequencies 
(Mueller et al. 2014, Roelfsema et al. 2012), and where A and 
(Tnib are both expressed in /jm. As shown in Table A. 1, //mb.o and 
(Tnib are different for the different mixers of HIFI. 


3. Results 


We detected emission from the molecular ions HCO^, 
H13C0+, SO+, N 2 H+, and (tentatively) H 3 O+ toward OH 231.8 
(Table 1 and Figs. 2-5). Except for HCO^ 7=1-0 (Morris et 
al. 1987; Sanchez Contreras et al. 1997), HCO^ 7=3-2, and 
H^^CO^ 7=1-0 (Sanchez Contreras et al. 1997), all transitions 
reported here are first detections in this source. SO^ and HsO^ 
have not been detected before in CSEs around evolved stars 
(either oxygen- or carbon-rich). Quintana-Lacaci et al. (2013) 
recently reported the discovery of N 2 H^ emission toward the 
massive yellow hypergiant IRC-h 10420, but no detection of 
this molecule in O-rich low-to-intermediate mass evolved stars 
(AGB or post-AGB) has been published to date. Detection of 
N 2 H^ has been reported in two C-rich objects, the young PN 
NGC7027 and the pre-PN CRL618 (see § 1). 

In Table 1 we list the main parameters of the molecular tran¬ 
sitions observed, which are all in the ground vibrational state. 
The flux is computed by integrating the area under the line pro¬ 
file. The corresponding formal errors of the flux are given in 
parenthesis in the last column of the table and have been calcu¬ 
lated according to the error propagation theory as 


Table 1. Parameters of the transitions used in this research: 1) rest fre¬ 
quency (MHz); 2) rotational quantum numbers; 3) upper level energy 
(Eu); 4) spontaneous emission Einstein’s coefficient (A^i); and 5) inte¬ 
grated intensity or line flux (f T^dv) and its error (within parenthesis); 
for nondetections, we provide 3cr-upper limits to the line flux adopt¬ 
ing a spectral window of 30km s“', except for for which we use 
12 km s“' (see § 3.4). The permanent dipole moment (p) of the ions are 
also indicated. 


Rest Freq. 
(MHz) 

Transition 

QNs 

E„ 

(K) 

Aui 

(s-‘) 

jr; 7 u 
(Kkm s ') 

H 

89188.5 

1-0 

3.888 Debyes 

4.3 4.161E-05 

4.56 (0.05) 

267557.6 

3-2 

25.7 

1.444E-03 

9.1 (0.30)* 

535061.6 

6-5 

89.9 

1.244E-02 

0.15'' (0.05) 

624208.4 

7-6 

119.8 

1.997E-02 

<0.11 

H 

86754.3 

' 3 CO+ p= 

1-0 

3.888 Debyes 

4.2 3.829E-05 

0.50 (0.05) 

260255.3 

3-2 

25.0 

1.329E-03 

1.85 (0.11) 

346998.3 

4-3 

41.6 

3.267E-03 

0.76 (0.10) 

520459.9 

6-5 

87.4 

1.145E-02 

<0.11 

162198.6 

SO^ ju- 2.3 Debyes 

7/2-5/2(e) 16.7 

1.108E-04 

0.64 (0.04) 

208965.4 

9/2-7/2(f) 

26.8 

2.457E-04 

0.50 (0.08) 

254977.9 

ll/2-9/2(e) 

38.9 

4.566E-04 

0.55 (0.15) 

301361.5 

13/2-ll/2(e) 

53.4 

7.655E-04 

<0.15 

301736.8 

13/2-ll/2(f) 

53.5 

7.684E-04 

<0.15 

347740.0 

15/2-13/2(e) 

70.1 

1.189E-03 

0.25+ (0.15) 

348115.2 

15/2-13/2(0 

70.2 

1.193E-03 

<0.17 

486837.2 

20/2-18/2(0 

133.5 

3.330E-03 

<0.11 

487212.1 

20/2-18/2(e) 

133.7 

3.338E-03 

<0.11 

93173.1 

N2H+ p= 

1-0 

3.4 Debyes 
4.5 

3.628E-05 

0.54 (0.04) 

279510.8 

3-2 

26.8 

1.259E-03 

1.17(0.10) 

558965.1 

6-5 

93.9 

1.085E-02 

<0.12 

652094.2 

7-6 

125.2 

1.741E-02 

<0.14 

307192.4 

H3O+ n= 

1--2| 

1.44 Debyes 
79.5 

3.498E-04 

0.28 (0.09) 

984708.7 

Oo - lo" 

54.6 

2.305E-02 

0.40+ (0.20) 

1031299.5 

4+ - 3" 

232.2 

5.148E-03 

<0.31 

1069827.6 

4 + . 3 - 

^7 ■^7. 

268.8 

9.850E-03 

<0.23 


Elux after removing the contribution of the SO 2 line (~50%, see 
§3.1). 

Tentative detection. 


not include additional uncertainties that may result from baseline 
subtraction and absolute flux calibration. Typically, we binned 
our original spectra to match a common velocity resolution of 
At) ~2km s ', except for the 2MHz-resolution WILMA data at 
frequencies lower than ~300 GHz, for which Au is already larger 
than 2 km s“' (§ 2 ). 

We obtained the line frequencies and other spectroscopic pa¬ 
rameters discussed here from the Cologne Database for Molecu¬ 
lar Spectroscopy (CDMS; Muller et al. 2005), the Jet Propulsion 
Laboratory (JPL) molecular spectroscopy database (Pickett et al. 
1998), and the MADEX catalogue (Cernicharo 2012). 


cr = rms VEWZI x At), 

where EWZI is the spectral window (in km s ') over which 
the line is integrated and rms is the root mean square statistical 
noise (in kelvins) of the spectrum near the line for a spectral 
resolution of Au (in km s '). This estimate of the flux error does 

^ The HIFI Beam: Release #1. Release Note for Astronomers at 
http://herschel.esac.esa.int/twiki/bin/view/Public/Hifi- 
CalibrationWeb#HIFI_performance_and_calibration 


3.1. H'^CO^ andH'^CO^ 

The three HCO^ lines we observed toward OH 231.8 are shown 
in Eig. 2. The single-dish profiles of the 7=1-0 and 7=3-2 tran¬ 
sitions were already known from previous observations (Sanchez 
Contreras et al. 1997,2000). As expected, we find similar shapes 
and absolute intensities (within the calibration uncertainties). 
The HCO^ 7=1-0 transition, observed over a full velocity range 
of ~ 200 kms“* (yLSR=[-80:H-120]kms“*), is characterized by 
a broad flattened profile centered around Vlsr^ 30km s * with 
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a full width at half maximum of FWHM~90kms ' and a 
weaker, blue-shifted component centered at Vlsr— 53 km s ' 
withFWHM~25kms-'. 

The HCO^ J-3-2 line is also broad, covering the velocity 
range VLSR=[-20:H-105]km s“*. There is a prominent feature 
centered at yLSR~55km s“' with FWHM~24km s“*, which rep¬ 
resents roughly half of the total flux measured in the profile. 
This feature, observed but not discussed in previous works, is 
a blend of HCO^ J-3-2 with the SO 2 /A:,yt=133,ii-132,i2 line at 
267.5 GHz (Eu=105K). SO 2 shows a wealth of intense transi¬ 
tions that dominate the mm/far-IR spectrum of OH 231.8. We 
have used these other SO 2 transitions to estimate the total col¬ 
umn density and excitation temperature of this molecule with the 
classic population diagram method (to be published in Sanchez 
Contreras et al., in prep). Using this technique, we confirm 
the identification of the SO 2 line as the main contributor to the 
FLSR~55kms“* feature and remove its emission from the to¬ 
tal flux measured. We find that the expected intensity of this 
SO 2 transition accounts for 50% of the total flux measured in the 
blend at 267 GHz (f r^t/r)~18 Kkm s ') and, thus, the remain¬ 
ing ~50% is due to HCO^ J-3-2 (Table 1). A two Gaussian 
fit to the blended profile is consistent with the HCO^ line be¬ 
ing centered at yLSR~40km s * and having FWHM~80km s * 
(although these figures are rather uncertain). 

Regarding HCO^ 7=6-5, although the low S/N ratio pre¬ 
vents an accurate characterization of the profile, the shape 
of this weak line is also consistent with a broad pro¬ 
file (FWZI~60km s *) with two peaks around yLSR~20 and 
50km s“*, similar to the profile of the H'^CO^ 7=4-3 (also in 
Fig. 2, right column). 

As already noted and discussed by, e.g., Morris et al. (1987) 
and Sanchez Contreras et al. (1997), the broad HCO^ profiles are 
notably different from those of most molecules, including '^CO 
and *^CO. The latter display an intense narrow emission core (in 
the range yLSR=[+10:-i-55]km s“') that arises at the low-velocity, 
low-latitude parts of the nebula (clump 13), which contain ~2/3 
of the total nebular mass, plus much weaker broad wings pro¬ 
duced by the less massive but faster bipolar lobes. (Profiles of 
'^CO transitions from our survey are included as additional ma¬ 
terial in AppendixB, Fig. B.2). 

The spatio-kinematic distribution of HCO^ 7=1-0 is known 
with ~4"x2"-resolution from interferometric maps by Sanchez 
Contreras et al. (2000) - see Fig. 1. These data show that, indeed, 
the different line profiles of HCO^ and CO come from significant 
differences between the spatial distribution of these molecules. 
In contrast to CO, the HCO^ 7=1-0 emission is enhanced in the 
bipolar lobes and it shows a deficit of emission from the massive 
slow central region; this explains the broad flattened HCO^ pro¬ 
files. The HCO^ emission peaks at a compact region (referred 
to as clump 12^) at ~3-4" along PA~50° from the center. There 
is a large velocity dispersion, of ~80km s"', within this clump 
that, indeed, produces the bulk of the HCO^ 7=1-0 emission in 
the velocity range yLSR=[-20:H-60]km s“'. The extent of the 
HCO^ emission in the outflow is smaller than that of CO (in the 
southern lobe, the HCO^ emission abruptly falls beyond clump 
15) although the axial velocity gradient is similar in both cases. 

The single-dish profile of the 7=1-0 transition (Fig. 2) is 
consistent with the whole HCO^ emission completely fill¬ 
ing the telescope beam at 89 GHz, HPBW=27'.'3: the bluest 
HCO^ 7=1-0 spectral feature at yLSR=-53km s"' traces the tip 
of the approaching North lobe (clump II), while the reddest wing 
emission (up to yLSR=+120km s“*) originates in the South lobe 
(mainly from clump 14 and partially from clump 15). 


The smaller full width of the HCO^ 7=3-2 transition 
(~120km s '), compared with the 7=1-0 line, is consistent with 
the smaller beam size of the IRAM 30 m telescope at 267 GHz 
(HPWB=9'.'l) and, thus, the smaller fraction of the axial outflow 
lying within the beam. Given the velocity range over which the 
HCO^ 7=3-2 emission is observed in our single-dish data, most 
of the emission must arise at the nebular center, the bright clump 
I2d and, partially, at the base of the South lobe. 

The full width of the HCO^ 7=6-5 line is ~60kms '. 
At this frequency, the beam of Herschel is sufficiently large 
(HPBW=39'.'8, i.e., larger than the IRAM 30 m beam at 89 GHz) 
as to include the bipolar HCO^-flow in its full extent. Therefore, 
one would expect to observe HCO^ 7=6-5 emission over a ve¬ 
locity range comparable to that of the 7=1-0 line (~200km s '). 
The narrower 7=6-5 profile could indicate that the emission 
from this higher excitation transition, which has Eu=89 K, is re¬ 
stricted to regions closer to the nebula center, with smaller ex¬ 
pansion velocities. In particular, the emission could mainly be 
associated with the bright emission clump Ej. However, these 
conclusions are very tentative because of the low S/N of this 
transition. 

The H'^CO^ line profiles are shown in Eig. 2 (right column). 
As expected, these transitions are weaker than their main iso¬ 
tope analogues. The spectra are broad, with a full width that 
ranges from ~190km s * in the 7=1-0 transition to ~70km s * 
in the 7=4-3. They lack a dominant narrow core component 
but, instead, show broad, multiple-peaked shapes. Eor the 7=1- 
0 transition, we clearly detect emisson at yLSR=-50km s“' from 
the tip of the North lobe (clump II). In all transitions, we detect 
emission from a velocity range that is consistent with emission 
produced at the bright clump 12^ and at the base of the South 
lobe (clump 14). The emission contribution of the latter, clump 
14, concentrates in the range Vlsr^[ 60:80] km s“*, this spectral 
feature being rather prominent and sharp in the 7=3-2 and 4-3 
lines. 

3.2. SO+ 

SO^ is a reactive radical with a ^H-electronic state. The energy 
levels of a ^H-molecule exhibit the additional splittings due to 
the electron spin and orbital angular momentum interactions (Q- 
and A-doublets). The rotational levels are defined with three 
quantum numbers, G, 7, 1, which represent the absolute value of 
the projection of the total electronic angular momentum on the 
molecular axis, the total angular momentum from rotation and 
electronic motion, and the parity, respectively (Amano, Amano, 
& Warner 1991). 

A total of twenty SO^ transitions (five A-doublets in each 
of the Q.-1/2 and Q=3/2 states) lie within the frequency range 
covered in our mm-wave survey with IRAM 30 m. The SO^ 
transitions detected toward OH 231.8 are listed in Table 1 and 
their profiles are shown in Eig. 3. The lines detected correspond 
to the lowest rotational levels of the ^ni /2 ladder of SO"^, which 
require low- excitation temperatures, between Eu=17 and 70 K. 
The nondetection of the ^n 3/2 transitions is in very good agree¬ 
ment with the expectations given that the lowest levels are at 
Eu^500 K above the ground, much higher than the average tem¬ 
perature (Tkin^lO-TOK) of the molecular outflow of OH231.8 
(§ 1 and § 4.1). In Table 1 we have not included additional SO^ 
transitions that may be tentatively detected but that appear par¬ 
tially or fully blended with other stronger lines or in particu¬ 
larly noisy spectral regions (the spectra over the full range of 
frequencies covering the five SO^ doublets observed are shown 
in Eig. B.l of the Appendix). 
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The SO^ line profiles are broad, with FWHM~60+10km s ' 
and centered around Vlsr^TO km s“* (Fig. 3). The moderate S/N 
of the spectra does not justify an accurate analysis of the profiles, 
however, it is sufficient to establish that, as for HCO^, the SO^ 
lines lack a prominent, narrow central core, but instead are dom¬ 
inated by a broad spectral component. This suggests a similar 
distribution for SO^ and HCO^ and, therefore, it is likely that 
the emission from both molecules arises mainly at the fast out¬ 
flow. 


3.3. N 2 H+ 

N 2 H^ is a linear, triatomic molecule. Its spectral lines show hy- 
perfine structure due to the electric quadrupole moment of the 
nitrogen nucleus. The line widths of the individual hyperfine 
components, which are separated by a few km s ', are expected 
to blend because of the large expansion velocities of OH 231.8’s 
molecular outflow. In Table 1 we list the N 2 H^ 7=1-0 and 3-2 
transitions detected in the mm-wavelength range and also pro¬ 
vide flux upper limits for nondetections from our far-IR survey 
with Herschel up to upper level energies of Eu~ 160 K. 

The spectra of the two N 2 H^ transitions detected are shown 
in Fig. 4. Both lines have rather broad, somewhat struc¬ 
tured profiles: the 7=1-0 line is centered at 45+4kms * 
and has FWHM=85+10km s * and the 7=3-2 line is centered 
at 39+4kms * and has FWHM=70+10km s *. The smaller 
line width of the N 2 H^ 7=3-2 transition is consistent with the 
smaller beam size at this frequency (HPBW~8'.'7), which im¬ 
plies that a smaller fraction of the outflow is probed. 

The deconvolved width of the lines at half maximum, taking 
the hyperfine structure and the spectral resolution of the data into 
account, is ~70km s“', indicating that most of the N 2 H^ emis¬ 
sion arises in the fast outflow. As for HCO^ and SO^, N 2 H^ 
lines present broad profiles, lacking of a prominent, narrow core 
component. 


coincides with the systemic velocity of the source within errors, 
which reinforces the detection case. 

We observed the o-H 30 ^ 7^=0q- 1| line at 985 GHz and 
other transitions of this molecular ion at higher frequencies 
in our survey with Herschel toward OH 231.8. The spectrum 
at 985 GHz (shown in Fig. 5) is consistent with weak emis¬ 
sion centered at yLSR=37+3km s“* and FWHM=8+5km s“*. 
The line flux measured over a 14km s '-wide window is 
0.4(+0.2)Kkm s“', so we take a conservative position and con¬ 
sider this as a tentative detection. Nevertheless, it is worth not¬ 
ing that very weak emission from this transition is indeed ex¬ 
pected if H 30 ^ arises in the warm (sslOOK) envelope compo¬ 
nent that we recently discovered in our Herschel survey (see 
§ 1 and also §4.1). This warm component is found to be se¬ 
lectively traced by certain molecules, in particular, by H 2 O (as 
well as by CS, H 2 S, and SiS). Since the chemistry of H 30 ^ and 
H 2 O are expected to be intimately related to each other, we con¬ 
sider the relatively high temperature suspected from the weak¬ 
ness of the H 30 ^ 985 GHz transition (compared to the H 30 ^ 
line at 307 GHz) consistent with both species being produced in 
similar regions. 

4. Analysis 

4.1. Column densities and excitation temperatures 

We have obtained beam-averaged column densities (Ntot) and ex¬ 
citation temperatures (Tex) for the different ions observed using 
the standard population diagram technique (e.g.. Goldsmith & 
Ganger 1999). In this method, the natural logarithm of the col¬ 
umn density per statistical weight (N„/g„) is plotted against the 
energy of the upper level above the ground state (E„) for a num¬ 
ber of transitions of the same molecule. Assuming that the lines 
are optically thin and thermalized, i.e., all levels are under local 
thermodynamic equilibrium (LTE) conditions at a given unique 
temperature, N„/gu and Eu are related by the following formula: 


3.4. H 3 O+ 

The molecular ion H 30 ^, hydronium, is isoelectronic with am¬ 
monia (NH 3 ) and has a similar physical structure: H 30 ^ is a 
trigonal pyramid with the oxygen atom at the apex and the hy¬ 
drogen atoms at the base. Like ammonia, H 30 ^ undergoes in¬ 
versions as the oxygen atom tunnels through the plane of the 
hydrogen atoms leading to the characteristic inversion splitting 
of its rotational levels. An energy level diagram of this molecule 
for the lowest four rotational levels is shown by, e.g., Mamon et 
al. (1987). The rotation-inversions transitions are described by 
three quantum numbers 7^, where the parity number p takes a -H 
or - value (see, e.g., Kozlov & Levshakov 2011). H 30 ^ has both 
ortho- and para- modifications: transitions with K-hn, where n 
is an integer > 0 , are from the ortho and the rest are para variety. 

The spectrum of the p-H 30 ^ 7^=1 j^-2| line at 307 GHz is 
shown in Eig. 5. A weak emission line is observed centered at 
FLSR=34+3km s“' and with a width of EWHM~14+6km s"'. 
This profile is significantly narrower that those of the other 
molecular ions described here, but consistent with the low ex¬ 
pansion velocities found in the central regions of OH 231.8 (§ 1). 
Although the H 30 ^ emission is weak, it satisfies a simple objec¬ 
tive detection criterion: more than three adjacent channels of 
width 3.9km s ', corresponding to an anticipated line-EWHM 
of ^12km s“', are above 1.5cr (this is adopting a gaussian pro¬ 
file with a ^3cr level at the peak). Moreover, the line centroid 


, Nu . 3kWei , Atot Eu 

In — = In —3-T = In-, 

gu Z(Tex) kTex 


( 3 ) 


where k is the Boltzmann constant; Wui is the source bright¬ 
ness temperature integrated over velocity; Vui and Sui are the fre¬ 
quency and line strength of the transition, respectively; p is the 
appropriate component of the permanent dipole moment of the 
molecule; Z(Tex) its partition function; and u and I refer to the 
upper and lower levels involved in the transitions. According to 
eq. 3, for a given molecule a straight-line fit to the points in the 
population diagram provides Atot/Z(7’ex) from the y-axis inter¬ 
cept and Tex from the slope of the fit. The line flux W is given 
by 


W-f,VeUndv, ( 4 ) 

where rjes is the ratio between the main-beam efficiency and the 
forward efficiency of the telescope at a given frequency (values 
are given in § 2) and /d is the dilution correction factor. The latter 
is estimated as 


/d = 1 - exp{- 


Gsrc 

Gmb 


Gsrc — 

Gmb = 


-M 


41n2 


HPBW'' 


where Qsrc and Gmb are the solid angular extent of the source and 
the telescope main beam, respectively. These solid angles are ex¬ 
pressed above as a function of the telescope beam (HPBW) and 
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the angular major ( 0 a) and minor ( 0 b) axes of a uniform elliptical 
source (see, e.g., Kramer 1997), which has been taken to roughly 
represent the elongated molecular outflow of OH 231.8. 

Based on the angular extent of the HCO^y=l-0 emission in 
OH 231.8 deduced from interferometric mapping (Sanchez Con¬ 
treras et al. 2000); see discussion in § 3.1, we adopt a characteris¬ 
tic angular size of 0a x 0b ~4"xl2" for all the species discussed 
here. Although the spatial distribution of H^^CO^, SO^, and 
N 2 H^ is unknown, their relatively broad emission profiles and 
similar excitation requirements justifies our assumption. In the 
case of H 30 ^, the narrower line profile and possibly higher Tex 
inferred (see below) may suggest a more compact distribution. 
We have also considered a smaller size of the H 30 ^-emitting 
area of 0a x 0b ~ 4"x4", which leads to a lower Tex but a com¬ 
parable value for Niot- We note that the adopted angular size 
affects both Nu/gu and on the one hand, the smaller the size, 
the larger the value of Nu/gu obtained from the fit; on the other 
hand, /d changes by a different factor for different transitions, 
which translates into a modification of the slope of the straight- 
line fit to the points in the population diagram and, thus, into a 
different In the particular case of H 30 ^, both effects com¬ 
pensate to yield a similar value of Ntot when adopting a smaller 
size of ~4"x4". 

In Fig. 6 we plot the population diagrams for the molecular 
ions detected in OH 231.8. The partition function, Z(7’ex), has 
been computed for each molecule by explicit summation of a 
finite number of rotational energy levels from the ground vibra¬ 
tional state using MADEX (Cernicharo 2012). For H 30 ^ we 
have adopted an ortho-to-para ratio of 1 , which is the expected 
value in statistical equilibrium conditions at T>5QK (Phillips et 
al. 1992). The values of Atot and Tex obtained for the differ¬ 
ent species are given in Table 2. In this table, we also include 
results for '^CO obtained in Velilla Prieto et al. (2014) using 
the same technique and assumptions. As discussed by these and 
other authors, transitions are optically thin, or only mod¬ 
erately opaque at the line peak toward the densest central parts, 
and are expected to be thermalized over the bulk of the nebula, 
with high average densities of ~ 10 ^- 10 ® cm“^, and even larger 
in low-latitude regions closer to the nebular nucleus. 

For HCO+, H13C0+, SO+, and N 2 H+, the values of Tex ob¬ 
tained are low, ~8-20 K, in agreement with previous estimates of 
the kinetic temperature in the outflow, which is ~ 10 - 20 K in the 
fast lobes and, somewhat larger, ~10-40K in the low-velocity, 
low-latitude component® (see § 1). The excitation temperature 
inferred for H 30 ^, sslOOK, is significantly larger than for the 
rest of the molecular ions (Fig. 6 ). This value of Tex is very 
uncertain, however, the upper limit to the flux ratio between 
the 985 GHz and the 307 GHz transitions indicates temperatures 
certainy larger than 40 K for the emitting region; otherwise, the 
985 GHz line should have been much stronger, and thus well de¬ 
tected, given the relatively low energy of the upper state level of 
this transition (Eu~55 K) compared to that of the 307 GHz line 
(Eu~80K). Moreover, it is suggestive that the excitation temper¬ 
ature guessed for H 30 ^ matches very well that deduced for H 2 O 
and its main isotopologues, H 2**0 and H 2 ^^ 0 , from our recent 
survey with Herschel (Sanchez Contreras et al. 2014). Consid¬ 
ering that the chemistry of H 30 ^ is expected to be closely linked 
to that of H 2 O, the obtained result is probably more than an un¬ 
lucky coincidence. 


® These temperatures are derived from low-y transitions and, thus, cor¬ 
respond to the external layers of the central component, deeper regions 
are expected to be progressively warmer closer to the central star. 


On top of that, the narrow profile of the H 30 ^ 7 ^=l/- 2 | 
line at 307 GHz suggests that the emission arises in the dense 
central parts of the nebula where the expansion velocities are 
smaller and the temperatures are presumably higher. Support for 
this comes from the fact that the narrow line width and lowish 
LSR velocity we find for H 30 ^ is, remarkably, also observed 
for the (7, K) — (2,2) and (3,3) lines of NH 3 by Menten & Al- 
colea (1995). According to these authors, these lines’ profiles are 
“characterized by a narrow line width of < 15 km s ' (FWHM) 
and a centroid velocity around 30km s '", which is in excellent 
agreement with the values we find for H 30 ^ (see Section 3.4). 
This gives, first, credence to the reality of our H 30 ^ detection. 
Secondly, the NH 3 (2,2) and (3,3) lines arise from levels that are 
at energies 64 and 123 K above the ground state, respectively, 
and from their ratio Menten & Alcolea calculate a rotation tem¬ 
perature of ss 50 K from which they conclude that “the narrow 
component is emitted from a relatively warm region”. By associ¬ 
ation, this supports our conjecture that the H 30 ^ emission arises 
from a warm environment. In contrast, the lower excitation NH 3 
( 1 , 1 ) line, which arises from the molecule’s inversion-split para- 
ground state, shows broader emission than the (2,2) and (3,3) 
lines over a 73 km s“* (FWZI) velocity range with a (higher) 
centroid of 41 km s“*. 

For the values of Tex and Atot implied by the population di¬ 
agrams, all the transitions detected are optically thin (t<0.1). 
The largest optical depth (at the line peak) is found for the 
HCO^ 7=3-2 transition, for which t ~ 0.20. In this case, we ap¬ 
plied a small opactity correction factor of result¬ 

ing in an insignificant effect on the Tex and Atot deduced from the 
fit. The optical depth of the HCO^ 7=1-0 and 3-2 lines is only 
expected to be larger than 1 for column densities ^4xl0'"' cm“^; 
for the rest of the molecular ions under discussion even larger 
column densities, ^ 10 *® cm“^, would be necessary to produce 
optically thick emission. 

In the previous work by Sanchez Contreras et al. (2000), 
a large opacity for the HCO^ 7=1-0 transition is also ruled 
out based on the HCO^/H*^CO'^(7=1-0) intensity ratio, which 
agrees well with the low isotopic ratio estimated from 

CO and other molecules in the fast outflow, where the emission 
is known to be optically thin (see also § 4.2). 

4.1.1. Non-LTE excitation effects 

When the local density of molecular hydrogen (hhj) is insuffi¬ 
cient to thermalize the transitions of a given molecule, depar¬ 
tures from a linear relation in the population diagram are ex¬ 
pected. We have performed a non-LTE excitation analysis of 
the HCO^ and N 2 H^ emission. We have considered radiative 
and collisional processes and solved the statistical equilibrium 
equations under the large velocity gradient (LVG) approxima¬ 
tion using the code MADEX (Cernicharo 2012, and references 
therein). We have adopted collisional coefficients for HCO^ and 
N 2 H^ from Flower (1999) and Daniel et al. (2005), respectively. 
For SO^ and H 30 ^ there are no collisional coefficients available 
in the literature, and thus non-LTE calculations cannot be done. 

The critical densities of the mm-wave transitions of HCO^, 
H'^CO^, and N 2 H^ observed by us are in the range nciit~ 10 ®- 
10® cm“^. These values are comparable to the mean densi¬ 
ties in the central core component (^ 10 ®- 10 ^ cm“^) and moder¬ 
ately larger that the typical densities of the high-velocity clumps 
where the bulk of the wing emission arises (nH 2 ~ 10 ® cm“^). 
Note that the lowest H 2 densities in the molecular outflow of 
OH231.8, nH 2 ~ 10 "'- 10 ^ cm“^ are only found at very large dis¬ 
tances from the star in the south lobe (clumps 15 and beyond; see 
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Alcolea et al. 2001), but these distant, most tenuous and rapidly 
expanding regions (at yLSR^ 120 km s“') do not contribute to 
the emission prohles of the molecular ions under discussion, as 
shown in Figs. 2-5. 

For our non-LTE excitation calculations, we have adopted a 
representative density of na^-lQ^ cm“^. We hnd that for both 
HCO^ and N 2 H^, only mild LTE deviations are expected at 
this density, which result in values of Tex and Niot (as deduced 
from the population diagram) slightly lower than the real Fkin 
and Niot- This means that the temperatures and abundances 
quoted in Table 2 could be lower limits. We have explored 
the range of input Tkin and Niot that would still be consistent 
with the observations if the characteristic density of the emit¬ 
ting regions is nH 2 ~ 10 ^ cm“^. Eor both HCO^ and N 2 H^, we 
hnd that the population diagram could also be reproduced with 
a temperature of Tkin^BO-TOK; for these temperatures, the to¬ 
tal column densities would be Atot(HCO^)~9xlO'^ cm“^ and 
A);ot(N 2 H^)~ 1 . 5 xlO'^ cm“^. Therefore, the column densities 
deduced under non-LTE and LTE conditions do not differ sig- 
nihcantly (less than about 15%). 

4.2. Fractional abundances 

We computed the fractional abundances of the ions relative 
to H 2 , X-Niot/N}i^ as X=X('^CO)xAtot/Mot(^^CO) and adopted 
X('^CO)=5xlO“^, as estimated by Morris et al. (1987). This 
value of the '^CO abundance is in the high end of the typical 
range for O-rich stars; in the case of OH 231.8, this is consistent 
with the particularly low isotopic ratio, ~5-10, measured 

in this object (and other O-rich CSEs; Sanchez Contreras et al. 
2000; Milam et al. 2009; Ramstedt & Olofsson 2014, and ref¬ 
erences therein). We cannot exclude a lower value of 2f('^CO), 
by a factor ~2-3, in which case the abundances derived for the 
rest of the species should be scaled down in the same proportion. 
The fractional abundances obtained from the rotational-diagram 
method are listed in Table 2. As shown in § 4.1.1, the true values 
of X(HCO^) and X(N 2 H^) could be slightly larger than those 
given in the table if mild non-LTE excitation effects are present 
(discrepancies ^15% are unlikely). 

The average fractional abundance of HCO^ in OH 231.8 we 
derived, 2f(HCO^) sslO^*, is in good agreement with previous 
estimates (Morris et al. 1987; Sanchez Contreras et al. 1997, 
2000). Small differences are expected due to data calibration 
uncertainties and the slightly different 2 f('^CO) abundance and 
excitation temperature adopted in different works. The average 
abundance of HCO^ in OH 231.8 is also within the range of val¬ 
ues observationally determined on a small sample of O-rich en¬ 
velopes, ~0.15-1.3x10“^ (Pulliam et al. 2011). The abundance 
of H'^CO^, and its comparison with that of HCO^, is consis¬ 
tent with a low '^C/'^C isotope ratio of ~7, in agreement with 
previous estimates for this object (see above). 

The abundances of the other ions we presented have not 
been previously reported in O-rich CSEs. The fractional abun¬ 
dance of SO^ in OH 231.8 implies an abundance relative to SO 
of X(SO^)/ 2 f(SO)ss 2 xlO“^; this is taking into account the frac¬ 
tional abundance of SO in OH 231.8 computed in other works, 
X(SO)~[l-3]xlO-*’ (Guilloteau et al. 1986; Morris et al. 1987; 
Sanchez Contreras et al. 2000, and this survey, Sanchez Con¬ 
treras et al., in prep.). Eor N 2 H^, our abundance estimate is 
consistent with the upper limit provided by Morris et al. (1987) 
based on the nondetection of the N 2 H^ 7=1-0 line by these au¬ 
thors ( 2 f(N 2 H^)< 1.6x10“*). We deduce an average HCO^ abun¬ 
dance relative to N 2 H+ of X(HCO+)/X(N 2 H+)~ 6 ; as we will 
see in § 5.1, the previous ratio may represent an upper limit to 


Table 2. Results from the rotational diagram analysis (§4.1) and molec¬ 
ular fractional abundances, X, i.e., relative to H 2 , derived (§4.2). We 
have adopted ^(‘^COj^SxlO-^ (Morris et al. 1987). 


Species 

0aX0b 

Niot 

Tex 

X 


nxn 

(cm“^) 

(K) 



4x12 

3.0x10®' 

20 

5.0x10“'’ 

h'2co+ 

4x12 

8.0x10'* 

14 

1.3x10“* 

H'*CO+ 

4x12 

1.2x10'* 

10 

2.0x10“® 

SO+ 

4x12 

2.2x10'* 

20 

3.7x10“® 

N2H+ 

4x12 

1.3x10'* 

8 

2.2x10“® 

H3O+ 

4x12 

4x4 

~4xl0'* 

-5x10'* 

=sl00 

-55 

-7x10“® 


the [C0]/[N2] proportion in OH 231.8. Eor H 30 ^, we obtain a 
fractional abundance of ss7xl0“®, although we caution that this 
value is particularly uncertain given the low S/N of the HsO^ 
spectra. Our estimate is consistent with the range of values 
measured in interstellar clouds, 2f(H3O^)ssl0“'®-10“® (Woot- 
ten et al. 1991; Phillips et al. 1992; Goicoechea & Cernicharo 
2001). In these interstellar regions, the fractional abundance of 
H 30 ^ is expected to relate to that of water by a factor ssl/1000- 
1/6000, based on gas-phase model predictions and observations. 
If this proportion also holds in OH 231.8, then one would in¬ 
fer an order-of-magnitude water content of 2 f(H 2 O)ssl 0 “^- 10 “®. 
This value is at the low edge of the water abundance range de¬ 
duced for a sample of six evolved O-rich stars, 2f(H2O)~[10“*- 
2 x 10 “®], with the majority of the stars having 2 f(H 2 O)~ 6 xl 0 “"^ 
(Maercker et al. 2008). 

We recall that the abundances given in Table 2 represent av¬ 
erage values within the emitting regions covered by the telescope 
beam, that is, the central core and (partially) the fast lobes. The 
large wing-to-core intensity ratio of the HCO^ prohles is known 
to rehect an abundance contrast of a factor ~3-4 between the fast 
lobes and the slow central clump 13, as determined by Sanchez 
Contreras et al. (1997) from their analysis by individual spectral 
and spatial components of the HCO^ 7=1-0 emission maps; see 
their Table 2. This result is also corroborated by the interfero¬ 
metric observations of this transition (Pig. 1), which show that, 
in contrast to '^CO, '*CO, and most molecules, HCO^ is found 
in abundance in the lobes of OH 231.8, but it is scarce at the 
nebula center (Sanchez Contreras et al. 2000). 

The broad prohles of H'*CO^, SO^, and N 2 H^ suggest anal¬ 
ogous lobe-to-core abundance enhancements for these ions. Un¬ 
fortunately, except for the HCO^ 7=1-0 transition, obtaining 
column densities separately for the different velocity ranges cor¬ 
responding to the slow core and the fast lobes is hampered by 
the weakness of the molecular lines. However, a straightforward 
argument in support of a lobe-to-core abundance ratio >1 for 
H'*C0^, SO^, and N 2 H^ (as well as, of course, for HCO"^) is 
provided by the large equivalent-width velocity (Few) of these 
lines; the equivalent-width velocity is given by the ratio of the 
line Hux, J T^dv, to the peak intensity, (see, e.g., Sanchez 
Contreras & Sahai 2012). As can be seen in Pig. 7, the Few of 
the emission lines by these species are systematically larger than 
those of '*CO. The '*CO prohles are dominated by the most 
massive component, that is, the slow central core (clump 13). 
The larger Few of the H'*C0^, SO^, and N 2 H^ prohles indi¬ 
cates that the contribution to the emitting prohle from the less 
massive fast lobes is larger for these ions than for '*CO, imply¬ 
ing that their fractional abundances in the fast lobes are higher 
than in the slow massive clump 13. This also indicates that the 
true abundances of HCO^, H'*CO^, SO^, andN 2 H^ in the lobes 
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Table 3. Upper limits to the column densities and fractional abundances 
of other ions in OH 231.8 adopting a 4" x 12"-sized outflow, Tcx=l5 K, 
andFWHM=40kms-‘. 


Ion 

N( cm-A 

X 

Ion 

At(cm-2) 

X 

OH+ 

< 3 xl 0 ‘^ 

< 5 x 10 -'* 

HOC+ 

<1x10*2 

<2x10-*** 

H2O+ 

<2xl0'^ 

< 3 x 10 -'* 

H2COH+ 

<1x10*2 

<2x10-'* 

CO+ 

<2x10*2 

< 3 x 10 -*° 

H2NCO+ 

< 5 x 10*2 

<8x10-*** 

HCS+ 

<8x10*2 

<1x10-'* 

HSCO+ 

< 7 x 10*2 

<1x10-** 

SH+ 

<2x10*2 

< 3 x 10 -'* 

HOCS+ 

< 4 x 10*2 

< 7 x 10 -“* 

CS+ 

< 5 x 10*2 

<8x10-'* 

HC3NH+ 

< 4 x 10*2 

< 7 x 10 -“* 

CH+ 

<1x10*2 

<2x10-'* 

HOCO+ 

< 7 x 10*2 

<1x10-“* 

HCNH+ 

<6x10*2 

<1x10-** 

H2CI+ 

< 4 x 10*2 

<8x10-“* 

NO+ 

< 4 x 10*2 

< 7 x 10 -'* 

NeH+ 

<1x10*2 

<2x10-'* 

CF+ 

<8x10*2 

<1x10-'* 

H2F+ 

<2x10*2 

< 3 x 10 -'* 

C3H+ 

<2x10*2 

< 3 x 10 -'* 

ArH+ 

< 5 x 10*2 

<8x10-“* 

OH- 

< 3 x 10*2 

< 5 x 10 -'* 

NCO- 

<2x10*2 

< 3 x 10 -'* 

C2H- 

<1x10*2 

<2x10-*° 

C4H- 

<2x10*2 

< 3 x 10 -“* 

CN- 

<1x10*2 

<2x10-'* 

C3N- 

<1x10*2 

<2x10-“* 

SH- 

< 3 x 10 *'* 

< 5 x 10 -** 





and in the slow central clump are higher and lower, respectively, 
than the beam-averaged values given in Table 2. 

In the case of H 30 ^, the narrow profile of the weak 
2| transition at 307 GHz suggests that the bulk of the emission 
is produced in the slow central clump 13. The Vew of this tran¬ 
sition is very uncertain but if it is confirmed to be lower than 
that of '^CO (by higher S/N observations), as suggested by our 
data (Fig.7), then one would infer a larger abundance of this ion 
in the slow core than in the fast lobes (i.e., a lobe-to-core abun¬ 
dance ratio <1). Moreover, since the higher Tex suspected for 
compared to that deduced for '^CO favors the origin of 
the H 30 ^ emission in envelope layers closer to the center, the 
average abundance in Table 2 could represent a lower limit to 
the true value in the central, compact emitting regions. 

4.2.1. Upper limits to nondetections 

We provide upper limits to the column densities and fractional 
abundances of other ions that are not detected in OH 231.8 and 
that may be of interest (Table 3). We have adopted common val¬ 
ues for the source emitting size, 6^ x 9^ ~4"xl2", rex=15K, 
and FWHM=40km s ' and LTE conditions. For the ions con¬ 
sidered, adopting a smaller size, larger excitation temperatures, 
and narrower profiles (as may be the case of H 30 ^) normally re¬ 
sults in slightly smaller values of N and X than those given in the 
Table. For the different ions, the upper limits to A^tot are derived 
from the rms noise around the transition that is expected to have 
the largest S/N in our data (adopting K); we have always 

checked that the upper limits obtained are also consistent with 
the nondetections of the rest of the transitions in the frequency 
range of our observations. 

5. Formation of HCO+, SO+, N2H+, and H3O+ 

5.1. Chemical model 

We performed chemical kinetics models to investigate the for¬ 
mation of the observed ions in OH 231.8. Our model is based on 
that by Agundez & Cernicharo (2006) and we have used it re¬ 
cently to study the new N-bearing species detected in OH 231.8 
by Velilla Prieto et al. (2014). The code has also been em¬ 
ployed to model the chemistry in different astrophysical environ¬ 
ments, including the prototypical C-rich star IRCh- 10216 (see, 
e.g., Agundez et al. 2007, 2008, 2010b, 2012; Cernicharo et 


al. 2010, 2013), and the O-rich yellow hypergiant IRCH-10420 
(Quintana-Lacaci et al. 2013). The chemical network in our code 
includes gas-phase reactions, cosmic rays, and photoreactions 
with interstellar UV photons; it does not incorporate reactions 
involving dust grains, X-rays or shocks. Our code does not ex¬ 
plicitly consider the chemistry of isotopologues, that is, does not 
include isotopic fractionation, such as selective photodissocia¬ 
tion or isotope-exchange reactions; as a consequence, predic¬ 
tions for the H'^CO^ abundance distribution are not made. 

As input we adopt the same physical model for the envelope 
used in Velilla Prieto et al. (2014), which consists of two main 
components: a) a dense, spherically symmetric wind expanding 
at low velocity, which represents the slow, central nebular com¬ 
ponent of OH 231.8 (clump 13); and b) a rectangular slab of gas 
(plane-parallel geometry) with characteristics similar to those of 
the walls of the hollow bipolar lobes. We modeled separately 
these two components, which are externally illuminated by the 
interstellar UV field. 

The main parameters of the simple physical model for 
OH 231.8 used as input in our chemical code are summarized 
in Table 4; additional details justifying the temperature, density, 
and velocity laws adopted are given in § 6.1 of Velilla Prieto et 
al. (2014). The region of the envelope modeled in case a goes 
fromRin~10'^ cm (~20R*) to its end, at Rout~7xl0*® cm (we re¬ 
fer to this region as the intermediate/outer envelope). The outer 
radius is observationally determined (Sanchez Contreras et al. 
1997; Alcolea et al. 2001) and the inner radius has been chosen 
to be well beyond the dust condensation radius so that the full 
expansion velocity of the gas has been reached. In these inter¬ 
mediate/outer regions of the envelope, the chemistry is driven by 
chemical kinetics. Another major input of our chemistry model 
is the set of initial abundances of the ’parent’ species. These are 
formed in deeper layers and are injected to the intermediate/outer 
envelope. The parent species, typical of O-rich environments, 
and the initial abundances adopted in our model are as in Velilla 
Prieto et al. (2014) and are reproduced again in Table 5. 

As in Velilla Prieto et al. (2014), the sources of ionization 
and dissociation adopted in our model are cosmic-rays and the 
interstellar ultraviolet radiation field. The cosmic ray ionization 
rate adopted is 1.2x10 *’s * (Dalgarno 2006). The intensity 
of the UV field assumed is the Draine field or Go=1.7 in units 
of the Habing field (Go=l.6x10“^ergs ' cm“^; Habing 1968; 
Draine & Salpeter 1978). We computed the dust optical extinc¬ 
tion, Ay, across the envelope layers adopting the standard con¬ 
version from H 2 column density V/f2=9.3xlO^''xAy cm“^, de¬ 
rived by Bohlin et al. (1978). The gas-to-dust mass ratio implicit 
in this conversion, ~100, is comparable within errors to the aver¬ 
age value measured in nearby O-rich AGB stars and, in particu¬ 
lar, in OH 231.8 (e.g., Knapp 1985; Kastner & Weintraub 1995; 
Sanchez Contreras et al. 1998). 

The slow, central component. We used our chemical ki¬ 
netics model first to investigate the formation of HCO^, SO^, 
N 2 H^, and H 30 ^ in an O-rich AGB CSE similar to the slow, 
core component of OH 231.8 (model a, Table 4). The fractional 
abundances predicted by the model as a function of the distance 
to the center are shown in Eig. 8. Eor HCO^, N 2 H^, and H 30 ^, 
we checked that there is good agreement between our predictions 
and those obtained from analogous photochemistry simulations 
by Mamon et al. (1987); these authors modeled the chemistry of 
these three ions for O-rich CSEs with standard properties and a 
set of different mass-loss rates, including Mq yr“'. We 

found no previous modeling attempts for SO^ in circumstellar 
envelopes in the literature. 
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Table 4. Parameters of the central Mira star and molecular envelope of 
OH 231.8 used for the chemical kinetics models (§ 5.1.) 


Parameter 

Value 

Reference 

Distance (d) 

1500 pc 

b 

Stellar radius (R.) 

4.4x10*^ cm 

g 

Stellar effective temperature (T.) 

2300 K 

g 

Stellar luminosity (L,) 

10 " Eg 

g 

Stellar mass (M,) 

IMo 

g 

model a: slow, core 

AGB CSE expansion velocity (Eexp) 

20 km s ' 

a,e,f,g,i 

AGB mass loss rate (M) 

10 ^" Mo yr-‘ 

e,a,f 

Kinetic temperature (Tkin) 


i 

model b: lobe walls 

Wall thickness 

1 " 

a 

H 2 number density (MH 2 ) 

10 ^ cm“^ 

a 

Kinetic temperature (rpi„) 

20 K 

a 


Notes, a: Alcolea et al. (2001), b: Choi et al. (2012), c: Cohen (1981), 
d: Kastner et al. (1992), e: Morris et al. (1987), f: Sanchez Contreras et 
al. (1997), g: Sanchez Contreras et al. (2002), h: Sanchez Contreras et 
al. (2004), i: Velilla Prieto et al. (2014). 

Table 5. Initial abundances relative to H 2 for representative elements 
and parent molecules used as input for the chemical kinetic models. 


Species 

Abundance 

Reference 

He 

0.17 

a 

H 2 O 

3.0x10-" 

b,TE 

CO 

3.0x10-" 

c,TE 

CO 2 

3.0x10-’ 

d 

NH 3 

4.0x10-® 

e 

N 2 

4.0x10-® 

TE 

HCN 

2 . 0 x 10 -’ 

f,g 

H 2 S 

7.0x10-* 

h 

SO 

9.3x10-’ 

f 

SiO 

1 . 0 x 10 -® 

i 

SiS 

2.7x10-’ 

j 


Notes, a: Asplund et al. (2009), b: Maercker et al. (2008), c: Teyssier 
et al. (2006), diTsuji et al. (1997), e: Menten & Alcolea (1995), 
PBujarrabal et al. (1994), g: Schdier et al. (2013), h:Ziurys et al. 
(2007), i: Gonzalez Delgado et al. (2003), j: Schdier et al. (2007) TE: 
From LTE chemical calculations by Velilla Prieto et al. (2014)). 


As the gas in the envelope expands, parent molecules start to 
be exposed to the interstellar UV radiation and photochemistry 
drives the formation of new species in the outer layers. Pene¬ 
tration of photons through the envelope is gradually blocked by 
dust and also by self-shielding of abundant gas species (mainly 
H 2 , H 2 O and CO). At the deepest layers of the intermediate/outer 
envelope inaccessible to the UV radiation, cosmic rays are major 
drivers of the chemistry. Cosmic rays ionize molecular hydrogen 
(Hp , which combines with neutral H 2 to form . The latter is 
a fundamental ion for breaking down parent molecules leading 
to the formation of new ions in the inner layers of the envelope. 

It is well known that the radical OH is an important product 
of the photodissociation of H 2 O in the external regions of O- 
rich CSEs. Our model predicts that OH reaches maximum abun¬ 
dance at ~5xl0'® cm, in good agreement with the characteristic 
radius of the OH maser toroidal shell observed at the center of 
OH 231.8 (Morris et al. 1982; Zijlstra et al. 2001). Our model 


indicates that, as the radical OH, molecular ions should be most 
efficiently formed in the outermost layers of the envelope, show¬ 
ing a shell-like distribution with maximum abundances at typical 
distances from the star of ~[5-7]xl0'®cm. Farther out, there is 
an abrupt fall of the ions’ abundances due to their destruction by 
dissociative recombination with electrons (e~) and to photodis¬ 
sociation of their main progenitors. 

According to our model, HCO^ forms in the outer layers of 
O-rich envelopes mainly as a photodissociation product of CO 
and H 2 O (in order of importance): 

HOC+ -H H 2 ^ HCO-" -I- H 2 
CO+ -H H 2 ^ HCO+ 4- H 
C+ -H H 2 O ^ HCO+ 4- H ’ 

where the ion is predominantly produced in the CO photoion¬ 
ization chain, and the molecular ions HOC"^ and CO^ are formed 
through 

C+ -H OH ^ CO+ -H H 
CO+ -H H 2 ^ HOC+ 4- H • 

An additional, but nondominant pathway to HCO^ produc¬ 
tion is by proton transfer from N 2 H^ to CO (N 2 H'^ 4- CO —> 
HCO"^ 4 - N 2 ). The model predicts that HCO^ reaches high¬ 
est abundances in a very thin outer shell with a characteris¬ 
tic radius of ~ 6 xl 0 '®cm. HCO^ has, together with N 2 H^, 
the most sharply peaked distribution amongst the ions mod¬ 
eled. The model peak abundance in the shell is V(HCO^)sslO“^, 
that is larger than the average value derived in OH 231.8 
from these observations. Deeper into the envelope, i.e., in the 
~[10'^-3xl0'®] cm region before the inner rim of the peak 
abundance shell, HCO^ is formed in much lower amounts 
(V(HCO^)<10 ") directly from CO and H^, which is a prod¬ 
uct of cosmic-ray ionization (H 3 4 - CO —> HCO"^ 4 - H 2 ). 

The SO^ shell-like abundance distribution predicted by the 
model is broader than that of HCO^: for SO^ the inner rim of the 
shell is located deeper in the envelope than for HCO^. Within the 
SO^ peak abundance shell, the fractional abundance varies be¬ 
tween V(SO^)~2xlO“'^ and ~[3-4]xl0“*, the maximum value 
being reached at a distance of ~5xl0'® cm. The SO^ abundance 
distribution shows a “hump” or secondary peak at ~3xl0'®cm 
where the abundance is V(SO^)~7xlO“®. The values of the SO^ 
model abundance inside the shell seem, in principle, compatible 
with the beam-averaged value deduced in OH 231.8 from our ob¬ 
servations. The main formation paths of SO'^ in the outer shell 
are 

S+ 4- OH ^ SO+ -H H 
SO 4- hv —> SO'^4-e“ 

C+-H SO 2 ^ SO+-H CO • 


The reaction of with OH is the dominant formation route 
near the peak abundance radius; the hump or secondary peak 
of the SO^ abundance distribution results from direct photoion¬ 
ization of SO by UV photons and the reaction of C^ with SO 2 
(Eq.7). In regions deeper into the envelope, cosmic rays are the 
fundamental source of direct ionization of SO leading to the pro¬ 
duction of SO^ (although in much lower amounts, <10 than 
within the shell) and free electrons. 

The spatial distribution of N 2 H^ predicted by the model is 
sharply peaked, similar to that of HCO^, but with a very low 
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peak abundance of X(N2H^)~[0.9-1.0]xl0“*' that is reached at 
~[ 5 - 6 ]xl 0 ^®cm. This maximum abundance is more than two 
orders of magnitude lower than the mean value measured in 
OH 231 . 8 . Based on our model, N2H^ is formed mainly via 

H+ + N2 ^ N2H+ + H2 
N+ + H2 ^ N2H+ + H 
HCO+ + N2 ^ N2H+ + CO ’ 

where the first reaction above plays a major role all the way 
throughout the different layers of the intermediate/outer enve¬ 
lope, whereas the second and third reactions become an impor¬ 
tant factor in N2H^ production only in the outermost parts near 
the peak abundance radius. In these outer envelope layers, the 
ion N2 is mainly produced by the proton transfer reaction of 
He^, which is produced by direct ionization of He by UV pho¬ 
tons, with N2. 

Regarding H30^, our model predicts efficient formation of 
this ion in the outer envelope mainly as a product of the photodis¬ 
sociation chain of H2O and CO. The model peak abundance, 
X(H3O^)~5 x 10“*, which is larger than the observed value, is 
reached at ~ 5 xl 0 '®cm. Near the peak abundance radius, the 
distribution of H30^ is relatively broad and similar in shape and 
values to that of SO^. In these outer regions, the main reactions 
involved in the production of H30^ are 

H2O+ -I- H2 ^ H3O+ -H H 

HCO+ -I- H2O ^ CO -H H3O+ ( 9 ) 

In contrast to the other ions, H30^ is rather abundant not 
only within the peak abundance shell but also deeper into 
the envelope, with a fractional abundance that varies from 
X(H3O+)~2 x 10-“ at 10'^ cm to X(H3O+)~4xl0-® at the inner 
rim of the shell (at ~ 2 xl 0 ^®cm). In these inner envelope re¬ 
gions, H30^ is by far the most abundant of the ions modeled and 
is produced by protonation of water (H2O-1- H3 —» H30^ -H H2) 
after cosmic-ray ionization of H2 leading to H^. 

In the outer layers of the envelope, dissociative recombina¬ 
tion with electrons is the major mode of destruction of all the 
ions, except for N2H^, which instead is most rapidly destroyed 
by reactions with H2O and CO. Dissociative recombination with 
electrons ceases to be the main destruction mechanism for all 
ions deeper into the envelope, where the electron abundance is 
low (<fewxl 0 “^). In these inner regions, molecular ions are 
quickly disassembled primarily by reactions with H2O and NH3 
(and, in the case of N2H^, also with CO), i.e., 

HCO+ -I- H2O ^ CO -H H3O+ 

N2H++ H2O ^ N2 + H3O+ 

SO+-H NH3 ^ SO H-NH+ 

H3O+ H- NH3 ^ NH+ H- H2O 

The lobes. We have also modeled the chemistry in the walls 
of the lobes (model b), which are on average more tenuous than 
the central regions and, thus, enable deep penetration of the ion¬ 
izing interstellar UV radiation. The main input parameters of the 
simple physical model adopted for the lobe walls (a rectangular 
gas slab externally illuminated by one side by the interstellar UV 
field) are given in Table 4 . Our chemical kinetics code computes 
the evolution with time of the molecular fractional abundances 
inside a series of individual gas cells at different depths into the 


lobe walls. An example is given in Fig. 8b for a representative 
cell midway between the inner and outer edge of the walls (i.e., at 
~ 0'.'5 through the lobe wall). The depth of the cell into the wall 
is expressed as a function of the dust optical extinction. Ay. For 
the cell whose abundances are represented in Fig. 8b, the optical 
extinction is Ay~l mag. 

The spatial distribution of the model fractional abundances 
across the lobe walls attained at ~ 800 yr, which is the dynami¬ 
cal age of the outflow (§ 1 ), is shown in Fig. 8c. As a result of 
photodissociation, the abundances of all molecular ions rapidly 
increase inward across the lobe walls, reaching maximum values 
near Ay~l mag; at larger depths, the ions’ abundances gently de¬ 
crease as destruction processes by reactions with H2O, NH3, and 
CO start prevailing over formation routes. The peak abundances 
predicted by the model are comparable to the observed beam- 
averaged values in the case of SO^ , and possibly in the case 
of H30^ . In the latter case, the model value is lower than the 
observational value but the model-data discrepancy, of a factor 
~ 2 , is probably within the total uncertainties in the abundance 
calculation. According to the model, these two ions are the only 
ones that would form in detectable amounts^ in the fast lobes 
of OH 231 . 8 . In fact, in clear disagreement with the observa¬ 
tions, our model predicts H30^ to be the second most abundant 
ion (after SO^) in the fast lobes of OH 231 . 8 . For HCO^ and 
N2H^, the model predicts fractional abundances two orders of 
magnitude lower than observed. We have checked that the model 
abundances for these ions never reach values comparable to the 
observed values for a reasonable range of outflow ages, ~ 500 - 
lOOOyr. 

For all the ions discussed in this research, the maximum 
abundances reached within the lobes (model b) are smaller than 
the maximum abundances reached in the slow, central compo¬ 
nent (model a). The largest core-to-lobe peak abundance ratio 
(~ 500 ) is found for HCO^. The ion N2H^ is predicted to be 
almost equally deficient in the lobes and in the central, peak 
abundance shell. SO^ and H30^ are intermediate cases, with 
core-to-lobe peak abundance ratios of ~4 and 13 , respectively. 

It may appear a priori somewhat surprising that HCO^, 
which is predicted by the model to be efficiently produced 
in the outer regions of the slow, central component (with a 
peak abundance of V(HCO^)sslO“^), is so scarce across the 
lobe walls, where the maximum abundance reached is only 
V(HCO^)^ 10 “^°. This is because the HCO^ formation and de¬ 
struction processes in the lobes are similar to those prevailing 
in the deepest regions of the central component, where HCO^ 
is indeed not very abundant because of the rapid destruction of 
this ion by reactions with the abundant H2O molecule. In the 
peak abundance shell of the central component, however, the bal¬ 
ance between photochemical formation processes and dissocia¬ 
tive recombination with results in an efficient net production 
ofHCO+. 

5.1.1. Initial abundance of H 2 O. 

We have investigated how the abundances of the molecular ions 
under discussion depend on the initial abundance of H2O, which 
is a prime parent molecule tied up to the chemistry of these and 
other species. We ran our code using the same input parameters 
as in the previous section but with a relative abundance of water 
of 2f(H2O)=3xl0“^ that is ten times lower (Fig. 9 ). For a given 
ion, the final role of H2O either as a net builder or a net destroyer, 
after the balance of formation and destruction reactions with wa- 


^ For the sensitivity limit reached by our observations. 
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ter, varies throughout the envelope and, of course, is different for 
the different ions considered. 

We find that, within the peak abundance shell in model a, a 
decrease of the initial abundance of H 2 O results in a decrease 
of the peak abundances for all ions except for N 2 H^. The peak 
abundances of HCO^ and SO^ are lowered by a factor ~l/3-l/4 
and the abundance of H 30 ^ (and also that of the radical OH) 
by a factor ~1/10. In contrast, the peak abundance of N 2 H^ in¬ 
creases by a factor ~2-3. Moreover, in the low H 2 O abundance 
case, N 2 H^ no longer presents a sharply peaked, shell-like dis¬ 
tribution, but rather its abundance progressively increases from 
the inner layers up to a radius of ~[4-5]xl0'®cm; beyond this 
point, the abundance falls abruptly as for the rest of the species. 
The markedly different reaction of N 2 H^ to a decline in the H 2 O 
abundance with respect to the other ions is because, in the case 
of N 2 H^, H 2 O provides the main mode of destruction of this ion 
but is not a required formation ingredient in any of the regions 
of the envelope (see §5.1). In contrast, in the outer shell, for 
HCO^, H 30 ^, and SO^, water acts only as a main production 
agent (directly or indirectly via H 20 ^ and OH). 

Deeper into the envelope, HCO"^ and N 2 H^ are the two ions 
whose abundances are most sensitive to X(H 20 ): our model pre¬ 
dicts HCO^ and N 2 H^ enhancements of one order of magnitude 
for a similar decrease of the water abundance. In contrast to 
N 2 H^, the response of X(HCO^) to a change in the water abun¬ 
dance is different in the outer peak abundance shell and in re¬ 
gions deeper into the envelope. This is because, as explained 
earlier, in the inner envelope regions the formation of HCO^ is 
principally the result of protonation of CO from H^, with H 2 O 
acting almost exclusively as a main destroyer. We also find that 
the model abundances of SO^ and H 30 ^ vary only marginally 
in the inner regions when the H 2 O abundance is lowered. For 
SO^ this result is readily expected since neither the main for¬ 
mation nor destruction processes in the inner envelope involve 
water in this case (as discussed in § 5.1). The marginal response 
of H 30 ^ to a decrease of X(H 20 ), which a priori seems less pre¬ 
dictable than for SO^, is a direct consequence of the enhanced 
production of HCO^ in the inner envelope when 3 f(H 20 ) is low¬ 
ered. In the low-water case, because of the 3f(HCO^) increase, 
the dominant formation route of H 30 ^ in the inner envelope is 
by HCO^ -I- H 2 O —> CO + H 30 ^ (instead of water protonation 
by Hp. Therefore, the lower H 2 O abundance is compensated by 
the higher abundance of HCO^, which translates to an unaffected 
H 30 ^ abundance. 

In the lobes (model b), lowering the initial H 2 O abundance 
to X(H 20 )= 3 x 10 “^ produces a different effect at different depths 
into the lobe walls for HCO^, SO^, and N 2 H^. In the outer 
edge of the walls, fully exposed to the ISM UV radiation field 
(Av-Omag), the fractional abundances of these ions remain al¬ 
most equal or slightly smaller than in the high water abundance 
case (X(H 2 O)= 3 x 10 “'*), whereas in the shadowed (Av~2.3 mag) 
inner rim of the lobe walls the predicted abundances are higher, 
with the largest enhancements (of a factor ~20-40) being for 
HCO^ and N 2 H^. In the case of H 30 ^, a lower water abundance 
results in a lower abundance of this ion in all regions throughout 
the lobe walls; at the shadowed inner edge of the lobe walls only, 
both the high and low water abundance models predict a similar 
value of the H 30 ^ abundance of 2 f(H 3 O^)~ 2 xl 0 “®. 

For all ions modeled, it appears that lowering the water abun¬ 
dance results in the maximum abundance zone moving deeper 
into the lobe walls. As seen above, another effect is an increase 
of A(HCO^) and A(N 2 H^) in the lobes, however, the model pre¬ 
dictions are still significantly lower (~2 orders of magnitude) 
than the abundances observationally deduced. Moreover, the 


core-to-lobe abundance contrast for HCO^ and SO^ in the low 
water abundance case (~60 and 20 , respectively) is still irrec¬ 
oncilable with the broad profiles observed for these species. As 
in the high water abundance case, N 2 H^ is equally deficient in 
the core and in the lobes. In the case of H 30 ^, the core-to-lobe 
abundance ratio diminishes to ~2-3 for a low H 2 O abundance. 

Finally, by pushing the water abundance to the low end of 
the range of observed values in O-rich stars, A(H 2 O)ssl 0 “^, the 
abundances predicted by our chemistry model in all cases, except 
maybe for SO^, become much smaller than the observed values 
both in the slow core component and in the fast lobes. 

5.1.2. Enhanced elemental abundance of nitrogen 

In this section we examine the behavior of the ions discussed 
here against an overabundance of the elemental nitrogen (up to 
a factor x40)* and whether such an N-enhancement can explain 
the large fractional abundance of N 2 H^ observed in OH 231.8. 
As in Velilla Prieto et al. (2014), we ran our chemistry code again 
modifying the initial abundances of those parent species that are 
most sensitive to the elemental abundance of nitrogen, namely, 
N 2 , NH 3 , and HCN: N 2 is increased by a factor 40 and NH 3 and 
HCN by a factor 7 relative to the values in Table 5. We find that 
the N 2 H^ peak abundance increases to ~ 2 xl 0 in the slow¬ 
ing expanding shell and to ~10“^° in the fast lobes. Therefore, 
despite the substantial N-enhancement, the model still predicts 
N 2 H^ to be underabundant (by one order of magnitude). More¬ 
over, as in the case of a standard N abundance, the model fails 
to reproduce the enrichment of N 2 H^ in the fast lobes relative 
to the slow core deduced from the observations. This was, in 
fact, expected since changing the elemental abundance of nitro¬ 
gen comparably affects both the fast and slow outflow compo¬ 
nents. 

The elemental nitrogen enhancement has only a minor effect 
on the other ions discussed in this research, except for H 30 ^ and 
SO^ but only in the deepest layers of the slow central compo¬ 
nent. In these inner regions, the abundances of H 30 ^ and SO^ 
decrease (by a factor ~6 and ~ 10 , respectively) as these ions are 
destroyed at larger rates through their interaction with more NH 3 
molecules (eq. 10 ). 

Given the set of chemical reactions involved in the produc¬ 
tion and destruction of HCO^ and N 2 H^, the abundance ratio 
between these two ions provides an upper limit to the CO to N 2 
abundance ratio. In OH 231.8, we deduce A(HCO^)/A(N 2 H ^)~6 
from the observations and, therefore, one could conclude that 
[C 0 ]/[N 2 ]^ 6 . This value, which is slightly smaller than, 
but comparable to, the [CO]/[N 2 ]= 7.5 ratio in O-rich sources 
(adopted in our model. Table 5), does not rule out some N en¬ 
richment but it also does not clearly point to a significant ni¬ 
trogen excess. Based on observations and modeling of the 
[N ii]/l/l6548,6584A emission (and other optical lines) from the 
shock-excited regions at the tips of the lobes of OH 231.8, Co¬ 
hen et al. (1985) suggest that N may be overabundant in these 
regions by a factor ~5 relative to the solar value. This factor 
is much smaller than that used here to test this effect (x40), and 
that has been proven to be unable to bring the model N 2 H^ abun¬ 
dance close enough to the observed value. 

We also recall that an enhancement of the elemental nitro¬ 
gen is not able to reproduce satisfactorily the observed abun¬ 
dances of HNCO, HNCS, HC 3 N, and NO recently discovered in 
OH 231.8 (Velilla Prieto et al. 2014): on the one hand, HNCO, 


* Elemental N-enrichment can result from hot bottom burning (HBB) 
processes in the nucleus of ^3 Mq stars. 
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HNCS, and HC 3 N would be underestimated and, on the other 
hand, NO (and maybe others, such as NH 3 ) would be sig¬ 
nificantly overestimated. Moreover, there are other N-bearing 
molecules that, in contrast to N 2 H^, are not particularly abundant 
in OH 231.8. One example is HCN: while the average abundance 
of this molecule in O-rich CSEs is relatively large X(HCN)sslO“^ 
(Bujarrabal et al. 1994; Schoier et al. 2013), and actually larger 
than predicted by standard chemistry models (e.g., Nercessian 
et al. 1989), the value found in OH 231.8 is X(HCN)~fewxlO-* 
(Sanchez Contreras et al. 1997, 2014). 

In any case, an enhancement of the elemental N abundance 
would at most be able to reduce the model-data discrepancies 
for N-bearing species, but it is unlikely to help improve the 
model predictions for all the other molecules that are not chem¬ 
ical daughter products of N 2 , such as the ions HCO^, SO^, and 
H 30 ^ reported in this work. 

6. Discussion and conclusions 

Based on single-dish observations with the IRAM 30 m and Her- 
schel telescopes, we have detected HCO^, H^^CO^, SO^, N 2 H^, 
and (tentatively) H 30 ^ in OH 231.8. The broad line profiles 
(with FWHM of up to ~90km s ') and low-excitation temper¬ 
ature (-10-20K) deduced for HCO+, H'3C0+, SO+, and N 2 H+, 
support the location of these ions preferentially in the base of the 
cold high-velocity lobes of OH 231.8. The fact that the narrow 
emission component from the massive central parts of the neb¬ 
ula (clump 13) is minimal in the profiles of these ions, clearly 
points to an abundance enhancement in the fast lobes relative to 
the slow core. Although uncertain, the much narrower profile 
(FWHM~14km s ') and higher excitation conditions (sslOOK) 
suspected for H 30 ^ suggest a different spatial distribution for 
this ion, which may be more concentrated toward the central 
parts of the slow, central core. 

Our chemical kinetics model predicts that HCO^, SO^, and 
H 30 ^ can form profusely (with peak abundances of sslO^*) in 
the external layers of the slow central regions of OH 231.8 and 
thus, most generally, in similar O-rich CSFs. A moderate pro¬ 
duction of H 30 ^ (with fractional abundances of up to - 10 ^^) is 
also expected in regions deeper into the envelope. The produc¬ 
tion of N 2 H^ is, however, rather inefficient, leading to very low 
fractional abundances of 2 f(N 2 H+)$ 10 “". 

In the fast lobes of OH 231.8, our standard chemistry model 
predicts relatively high abundances for SO^ and H 30 ^ (of up 
to sslO^^); however, we predict that HCO^ and N 2 H^ are very 
scarce, with abundances of ^lO^'" and ;$ 10 “'^, respectively. 

There are various observational results that our chemistry 
model fails to reproduce satisfactorily. In particular, the high 
HCO^ and N 2 H^ abundances in the fast lobes deduced from the 
observed profiles. If these ions were present in the lobes in the 
low amounts predicted by the model, their emission would have 
remained undetected in our survey. Another model prediction 
that is inconsistent with the observations is that all ions should 
be more abundant in the slow core than in the fast lobes, how¬ 
ever, the broad profiles of HCO^, N 2 H^, and SO^, reflect exactly 
the opposite. One of the most extreme cases is HCO^, for which 
the model predicts a core-to-lobe abundance ratio of -500 versus 
the measured value of -1/3-1/4; the latter value is deduced from 
the broad, flattened profiles we report, and from accurate char¬ 
acterization of the HCO^ spatio-kinematic distribution obtained 
from previous single-dish and interferometric mapping (§ 1 ). 

Regarding H 30 ^, the model predicts a spatial distribution 
analogous to those of the other ions, i.e., peaking at the ex¬ 
ternal layers of the slow central core due to photodissociation 


processes. This distribution is not easily reconcilable with the 
suspected origin for the H 30 ^ emission in the deep (slow and 
warm) regions of the central core. The model predicts a de¬ 
tectable amount of H 30 ^, of ssl 0 “®, in the fast lobes that is not 
consistent with our observations. The tentative H 30 ^ abundance 
distribution in OH 231.8 deduced from our observations needs, 
however, to be confirmed by higher quality data. 

The large core-to-lobe abundance ratio of HCO^, SO^, and 
N 2 H^ predicted by the model can be reduced (but not sup¬ 
pressed) by adopting a lower initial H 2 O abundance of 3x10“^ 
(instead of 3x10 used in our standard model). Also, lower¬ 
ing the water content injected to the intermediate/outer envelope 
would produce somewhat lower values of the peak abundances 
of H 30 ^, which is probably overabundant in the external enve¬ 
lope layers in our standard model. 

We also find that an elemental nitrogen overabundance in 
OH 231.8 (by a factor of up to -40) is not able to produce suffi¬ 
cient N 2 H^, which remains one order of magnitude less abun¬ 
dant than observed. We believe that the model’s inability to 
explain the large abundance of N 2 H^ is not mainly a problem 
caused by an assumed nitrogen abundance that is too low. A 
similar problem also resides in NH 3 in the majority of O-rich 
CSFs, which contain ammonia in amounts that exceed predic¬ 
tions from conventional chemical models by many orders of 
magnitude (e.g., Menten et al. 2010; Menten & Alcolea 1995). 

6.1. Molecule reformation after the passage of fast 
dissociative shocks? 

Taking into consideration the various important model-data dis¬ 
crepancies deduced from our analysis, we must conclude that 
the standard molecular formation scenario adopted (triggered by 
ionization by UV photons and cosmic rays of gas whose ini¬ 
tial state is predominantly molecular) is too simple and/or does 
not correctly depict the recent molecular formation history of 
OH 231.8. The adopted scenario, however, may represent the 
chemistry of most molecules in the intermediate/outer layers of 
normal AGB CSFs reasonably well (see, e.g., Nejad & Millar 
1988; Glassgold 1996; Millar et al. 2000; Agundez et al. 2008). 

The notable chemical differences between OH 231.8 and 
most O-rich AGB CSFs are widely recognized and are not re¬ 
stricted to the ions discussed here but affect most of the species 
identified in this object (§ 1). Since the earliest studies on the 
molecular composition of OH 231.8, the remarkable molecu¬ 
lar richness of this source has been related to shocks. In¬ 
deed, the main difference between OH 231.8 and ’normal’ O- 
rich AGB CSFs is the presence of a fast (-400km s ') outflow 
in OH 231.8. If the acceleration of the bipolar lobes resulted 
from the violent collision between underlying fast jets on the 
pre-existing slow AGB wind, as proposed by many authors, then 
shocks must have necessarily played a major role in defining the 
current chemical composition of the outflow. 

In particular, given the high velocities observed in the fast 
outflow of OH 231.8, it is unlikely that molecules survived this 
kind of acceleration process but rather they were probably disso¬ 
ciated. Molecules are indeed expected to be destroyed by shocks 
with speeds larger than ^50km s ' (e.g., Hollenbach & McKee 
1980). After the shocked gas has cooled below -10"'K, atoms 
begin to associate and form molecules again; if the post-shock 
gas is moderately dense (e.g., sslO^- 10 ^ cm“^), molecule ref¬ 
ormation can happen indeed rather quickly, within a few years 
from the shock (e.g., Glassgold et al. 1989; Neufeld & Dalgarno 
1989). We believe that most molecules observed in the fast out¬ 
flow of OH 231.8, if shock accelerated, must have dissociated 
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and reformed in the post-shock gas. A similar scenario has been 
proposed to explain the chemistry in the high-velocity molecular 
outflow of the C-rich PPN CRL618 (Neri et al. 1992). 

At present, the shocks that accelerated the outflow ~800 yr 
ago are not active, therefore, the bulk of the molecular outflow 
in OH 231.8 is probably swept-up AGB wind material that has 
cooled down sufficiently to allow molecule reformation. Note 
that the lobe axial acceleration took place in less than sslSOyr, 
which may have been the duration of, e.g., a pulsed fast wind. 

The physical properties (density and temperature) of the fast 
post-shock gas, which could have been fully dissociated and 
fully ionized ~800yr ago, must have evolved as a function of 
time. In particular, the temperature in the lobes has certainly de¬ 
creased to its current value of ~10-20K. The assembly of new 
molecules would have happened progressively while the shock- 
heated material was cooling down. The ionized, atomic, and 
molecular fractions in the post-shock gas must have also evolved 
with time in a way that is not easy to predict, over the course of 
the cooling and molecular regeneration. [In particular, a transi¬ 
tion from an H^-dominated regime to an H 2 -dominated regime, 
passing through an intermediate H-dominated phase, has proba¬ 
bly happened.] Since these are the basic ingredients for putting 
together new molecules, their evolution with time must have dic¬ 
tated the set of chemical reactions that were dominant at each 
precise moment over the course of molecular regeneration lead¬ 
ing to the final molecular content observed at present. 

We increase the number of observational studies and chem¬ 
istry modeling of molecular ions, which is needed to make 
progress in determining the role of ion-molecule reactions in the 
chemistry scheme of O-rich CSEs. Taking the presence of fast 
outflows and the probable shock-acceleration history of the en¬ 
velope of OH 231.8 into account, the observations reported here 
may contribute to the emergence of an observational picture of 
the effects of shocks on the molecular chemistry of circumstel- 
lar envelopes. In the future, it would be desirable to attempt 
new chemical kinetics models adopting the molecule reforma¬ 
tion scenario proposed here for OH 231.8. This is a difficult task 
since determining the time dependence of the relevant physical 
parameters of the cooling gas, where second-generation species 
have formed, is a challenging theoretical problem. Undoubtedly, 
it would also be advantageous to enlarge the reaction network 
used in most standard chemistry models by including additional 
chemical reactions that may affect the predicted abundances of 
certain species (e.g., grain-surface reactions. X-rays, etc). 
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Fig. 7. Plot of the equivalent-width velocity (Vew) versus upper- 
level energy (E„) for HCO+, H‘-^CO+, SO+, N 2 H+, and H 3 O+ transitions 
(filled symbols) and *-^CO, for comparison (empty circles; profiles 
and line parameters are given in the Appendix). The dashed line is a fit 
to the Vew of the *^CO lines. For all ions, except for HsO^, the values of 
Vew are systematically larger than for '^CO, indicating a larger wing-to- 
core emission contribution to the observed profile; the H 30 ^ emission, 
in contrast, is dominated by the narrow, core emission component as 
indicated by the small Vew- 
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Fig. 2. Spectra of HCO"^ (left) and H'^CO^ (right) emission lines in OH 231.8. The velocity resolution (AV) and half power beam width (HPBW) 
are indicated inside the boxes at the top-right comer. A two Gaussian fit to the blend of the HCO^ 7=3-2 and SO 2 lines is shown (dotted lines). 
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Fig. 3. Same as in Fig. 2 except for SO^. See also figure B. 1 in the appendix. 
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Fig. 4. Same as in Fig. 2 except for N 2 H^. 
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Fig. 5. As in Fig. 2 except for H 30 ^. Detection of the o-H 30 ^line at 985 GHz is tentative. 
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Fig. 6 . Population diagrams. Line fluxes (circles) and upper limits (triangles) are from Table 1. The error bars include the rms of the spectra 
and absolute flux calibration uncertainties, up to ~25%. The frequencies of the transitions used are indicated in GHz. The column density and 
excitation temperature obtained from the linear fit to the data (dashed line) are given in Table 2. Article number, page 19 of 24 










Fig. 8. Results from our chemical kinetics model adopting a fractional water abundance of X(H20)=3xlO (§ 5 . 1 ). a) Spatial distribution of 

the model fractional abundances in the slow central core of the envelope (fossil remnant of the AGB CSE); b) evolution with time of the model 
fractional abundances in a representative gas slab or cell within the lobe walls (approximately in the midpoint between the outer and inner end of 
the wall) where the extinction is Ay ~1 mag; c) model fractional abundances as a function of the depth (expressed as Ay) into the lobe walls, at 
~800yr. Dotted horizontal lines in panel c) represent the observed beam-averaged values of our estimated fractional abundances. 
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Fig. 9. Same as in Fig. 8 except for a lower abundance of water, A'(H20)=3x 10 ^ - see § 5.1.1. 
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Appendix A: The HIFI beam pattern 

The HIFI beam pattern (size, efficiency, etc) is described by 
Roelfsema et al. (2012). In October 2014, a release note was 
issued by the HIFI calibration team (Mueller et al. 2014®) report¬ 
ing improved measures of the half-power beam width (HPBW) 
and main-beam efficiency (?7mb) that supersede previous esti¬ 
mates in Roelfsema et al. (2012). The new efficiency esti¬ 
mates are systematically lower than previous values by typically 
15-20% for T/nib- 

For convenience, we reproduce again in this appendix the 
updated values of the relevant beam-model parameters provided 
by Mueller et al. (2014) in Tables 2 and 3 of their release note. 
These values, given in Table A.l, are used to obtain the formulae 
(Eq. 2 in § 2.2) that describe the frequency/wavelength depen¬ 
dence of HPBW and /yeff for the HIFI bands (1 through 5) used 
in this work. Since we average the H and V polarization data to 
obtain final spectra, we use tfie average of /ymb for H and V at 
a given frequency in order to convert antenna temperatures (T'^) 
to a main-beam scale (Tmb) at that frequency. We recall that, as 
stated in § 2: 


Tmb - /Wff - 


m 

^mb 


Appendix B: Additional material 

In this Appendix, we include additional material for the ion SO^ 
and the molecule in OH 231.8. 

In Fig. B.l, we show all the SO^ transitions (detected and 
nondetected) observed within the frequency range covered in 
our survey with the IRAM 30 m telescope. Some of the SO^ 
fines are blended with transitions of other molecules and/or are 
in particularly noisy regions and have not been considered for 
the rotational diagram analysis presented in Section 4.1. 

A total of nine '^CO transitions have been detected in 
OH 231.8 from our surveys with IRAM 30 m and Herschel 
(Fig. B.2 and Table B.l). The ^^CO fines, wfiich are good tracers 
of the density and kinematics of the nebula, show structured pro¬ 
files with two main components: (1) a prominent, relatively nar¬ 
row core (centered around the systemic velocity of the source, 
y^j’^~33-34km s^b plus weaker, broad wings. The line core 
arises at the slow, dense central (low-latitude) parts of the neb¬ 
ula, which expand at relatively low velocities of ~10-35km s ' 
(clump 13), whereas the wing emission mainly arises at the fast, 
bipolar lobes (see § 1 and Fig. 1). (As explained in the Introduc¬ 
tion, the spatio-kinematic structure of the molecular outflow of 
OH 231.8 is well known from previous single-dish and interfer- 
emtric maps of '^CO). The large core-to-wing intensity ratio of 
the *^CO profiles is consistent with most (~70%) of the mass in 
the molecular outflow being in the slow central regions (Sanchez 
Contreras et al. 1997; Alcolea et al. 2001). 

The '^CO wing emission progressively weakens, relative to 
the core, at higher frequencies. For transitions observed with 
the same telescope (i.e., for the three 7u^3 fines observed with 
IRAM 30 m or for the six Ja>^ lines observed with Herschel), 
this trend is partially explained by the smaller beam and, thus, 
the smaller fraction of the fast outflow sampled by the observa¬ 
tions at higher frequencies. However, this is not the only rea¬ 
son, since at all frequencies the Herschel beam is comparable to, 

® The HIFI Beam: Release #1. Release Note for Astronomers - 
http://herschel.esac.esa.int/twiki/hin/view/Puhlic/Hifi- 
CalibrationWeh#HIFI_performance_and_calihration. 


Table A.l. Adopted values for HPBW, T/mi,, and Ruze-like scaling fac¬ 
tors T/mb.o and (Tmb (see Eq. 2 in § 2.2) for one spot frequency per mixer 
and per polarizations (H and V) of the HIFI instrument. Adapted from 
tables 2 and 3 of Mueller et al. (2014). 


Mixer 

/ 

(GHz) 

HPBW 

(") 

^mb 

^mb,0 

^mb 

ijim) 

IH 

480 

43.1 

0.62 

0.624 ± 0.007 

2.278 ± 0.004 

IV 

480 

43.5 

0.62 

0.618 ± 0.007 

2.248 ± 0.004 

2H 

640 

32.9 

0.64 

0.643 ± 0.009 

2.219 ± 0.003 

2V 

640 

32.8 

0.66 

0.662 ± 0.009 

2.224 ± 0.002 

3H 

800 

26.3 

0.62 

0.626 ± 0.008 

2.227 ± 0.004 

3V 

800 

25.8 

0.63 

0.639 ± 0.008 

2.287 ± 0.010 

4H 

960 

21.9 

0.63 

0.639 ± 0.008 

2.227 ± 0.004 

4V 

960 

21.7 

0.64 

0.643 ± 0.008 

2.252 ± 0.006 

5H 

1120 

19.6 

0.59 

0.595 ± 0.006 

2.068 ± 0.005 

5V 

1120 

19.4 

0.59 

0.597 ± 0.006 

2.114 ±0.005 


Table B.l. Same as Table 1 except for '^CO lines observed with IRAM 
30 m and Herschel. 


Rest Freq. 
(MHz) 

Transition 

QNs 

E„ 

(K) 

Aui 

(v') 

J 

(Kkm s“') 

1 

II 

o 

u 

0.11 Debyes 


110201.4 

1-0 

5.3 

6.336E-08 

13.45 (0.06) 

220398.7 

2-1 

15.9 

6.082E-07 

51.50 (0.13) 

330588.0 

3-2 

31.7 

2.199E-06 

42.00 (0.30) 

550926.3 

5-4 

79.3 

1.079E-05 

3.33 (0.07) 

661067.3 

6-5 

111.1 

1.894E-05 

3.03 (0.09) 

771184.1 

7-6 

148.1 

3.040E-05 

2.30 (0.13) 

881272.8 

8-7 

190.4 

4.574E-05 

2.00 (0.15) 

991329.3 

9-8 

237.9 

6.554E-05 

0.65 (0.15) 

1101349.6 

10-9 

290.8 

9.034E-05 

0.70 (0.15) 


or larger than, the IRAM 30 m beam, with a maximum size of 
HPBW=22'.'l at 110 GHz for the observations reported here. 

Both the full width of the wings and the full width 
at half maximum of the *^CO profiles decrease as the up¬ 
per energy level increases. The full width of the wings is 
largest for the *^CO (7=1-0) transition, which is observed 
over a velocity range of yLSR=[-80:-i-140]km s“* with a beam 
HPBW=22.1", and decreases for higher-7 transitions down to 
Flsr=[+15:-i-55] km s“* for the 7=7-6 line, observed with a 
comparable beam of HPBW=27.6"; for transitions with 7>7, 
we do not detect wing emission. The FWHM of the *^CO lines 
ranges between ~36 to ~14km s"' for the 7=1-0 and the 10- 
9 transitions, respectively. An analogous behavior is observed 
in most molecules observed in our surveys with IRAM 30 m 
and Herschel, including *^CO (Sanchez Contreras et al. 2014, 
Sanchez Conteras et al., in prep.). The observed trend sug¬ 
gests that the envelope layers with higher excitation conditions 
(i.e., warmer and, thus, presumably closer to the central star) are 
characterized by lower expansion velocities. 
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Fig. B.2. '^CO transitions detected in OH 231.8+4.2 with IRAM 30 m and Herschel. The red broad shoulder observed in the '^CO (7=8-7) line 
is probably an artificial feature due to residual baseline distorsion. 
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Fig. B.l. IRAM 30 m spectra toward OH 231.8 near the five mm-wave 
doublets of SO^ (see Table 1). An LTE model compatible with values 
in Table 2 and adopting an average FWHM=55km s“* for the lines is 
overplot (red line). 
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